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ABSTRACT 


Simulation of inventory systems is rapidly becoming an accepted 
management tool for quantifying proposed decision rules, The Stanford 
Research Institute's Inventory Management Simulator /Evaluator is an 
excellent example of a multi-branch, multi-item simulator, It represents 
one of the most sophisticated and comprehensive programs available today, 
In this Guide, the simulator is completely documented in such a way that 
anyone with a general knowledge of supply systems can readily follow 
the rationale of the program, For those who have a knowledge of the 
FORTRAN language and of programming techniques, a section is devoted 
to a card-by-card description of the program, Operating instructions 
and output reports are also discussed in detail. As an added feature, 
new sub-routines have been designed to produce random demand patterns 
which approximate the Normal, Exponential, Uniform, and Poisson distri- 


butions, 
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PREFACE 


My introduction to the Stanford Research Institute's Inventory 
Management Simulator/Evaluator occurred during a six week field assigne- 
ment to that organization in the summer of 1963, At the time I reported, 
Mr, Peter H. Butterfield, Manager of the Industrial Research Division, 
suggested that I work with Mr. Michael E. Chambreau on the redesign of 
their inventory simulator, Being interested in the inventory management 
field and fresh from a year of Operations Analysis studies I gladly 
accepted Mr, Butterfield's offer. The experience has proven very valuable 
and I am sure that it will be even more so in the years ahead, In par= 
ticular, the work of that summer helped to fill that gap which always 
exists between the theoretical knowledge of the classroom and the day=to= 
day application of that knowledge in the "real" world, It certainly gave 
me a better understanding and appreciation for both the theoretical and 
practical aspects of the operations researcher's work, Learning the art 
of applying a science is somewhat akin to putting the flesh on a skele= 
ton, You begin to recognize the value of the theoretical frame-work 
upon which a problem is solved as well as the value of the solution to 
that problem. It is with the general thought of helping others to 
bridge this gap that I have written this paper and highly recommend to 
succeeding students the study of actual research reports such as these 
listed in the bibliography. 

It is my hope that this volume will serve at least three specific 
purposes: first, to provide a publication which will give a detailed 


description of the computer program and necessary operating instructions 
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such that the simulator may be utilized by the students of the Operations 
Research, the Management, and the Management Data Processing Curriculums 
in their studies and thesis work; secondly, as a basic document on which 
to build a machine simulation course for Supply Corps Officers of the 
Operations Research curriculum which could be substituted for the present 
War Gaming course now required for these officers; thirdly, to provide a 
program which may be utilized as a model and source of basic building 
blocks with which to design more specialized programs, 

I would like to express my thanks to Mr. Butterfield for the courtesies 
he extended me during my visit at the Institute and for permission to use 
the simulator here at the Naval Fostgraduate School. Also, I am indebted 
to Mixe Chambreau for his patience and instruction in programming tech= 
niques and to Professor Rex Shudde and Commander S. W. Blandin for their 


advice and assistance on the preparation of this thesis. 
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CHAPTER JI 


INTRODUCTION 


In documenting the Stanford Research Institute's Inventory Manage= 
ment Simulator I have intentionally directed my efforts towards providing 
a working reference for students of the Naval Postgraduate School, An 
attempt has been made to explicitly point out the assumptions and ground 
rules of the program in order to define the system being simulated, It 
is felt that by clearly defining the system, tne simulator may be used 
by students who have no training in computer language as well as by those 
who have. I have also attempted to minimize the time required for a 
student to become thoroughly familiar with the program by presenting as 
many details as possible without overly boring the reader, The mathemati- 
cal derivations of the optimal requisition and reorder level equations 
have not been reproduced in this paper, Students who are interested in 
reviewing these are referred to reports 1 and 2 of the Bibliography. [I 
highly recommend to the Supply Corps Officers in the Operations Research 
curriculum the study of SRI's Technical Memorandum No, 3. In addition 
to providing the basis of the simulator program, it is also an excellent 
example of Operations Research at work, Chapter III which is a card-by= 
card description of the FORTRAN Language coding and program logic may be 


omitted by the student who is not familiar with the FORTRAN. 
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CHAPTER II 
THE SUPPLY SYSTEM TO BE SIMULATED 
l. General description 

The Supply System for which the simulator is designed, consists of 
a manufacturing facility, a central or "wholesale" warehcuse and a number 
of branch or "retail" warehouses, The system is intended to be a "push" 
system with a control office located at the central warehouse which 
periodically reviews at fixed intervals the system inventory positions 
and initiates action to resupply the system by placing production orders 
or resupplies a particular location by making interim shipments from ths 
central warehouse stocks, There is no provision made to reiocate excess 
stock at one branch to another which has a deficiency. 

The physical distribution of stock to the branches is made from the 
central warehouse which is located at the manufacturing center, Production 
orders are coordinated with the expected needs of the branches and alloca] 
tions are made to the branches at the time of completion of the manufac— 
turing process, The production quantity is physically delivered to the 
central warehouse which immediately computes the need of each branch based 
on its current inventory position and makes direct aliocations as required, 
Shipments are made to each branch on a fixed schedule depending upon the 
particular requirements of a branch, There are no provisions made for 
emergency type deliveries, The central warehouse acts as a branch or 
"retail" outlet within its region as well as being a stock point for ths 
entire system, 

The simulator was originally designed in connection with studies of 
products which have a common production operation up to a certain point 


and then are either modified slightly or packaged in different sizes which 





then become distinct items of the product in the inventory, These items 
retain a certain degree of interchangeability with respect to management 
decisions, There is, however, no provision in the simulator to substitute 
one item of stock for another to fill customer demands, However in re= 
viewing the system as a whole it is desirable to consider not only ths 
items individually but also to consider the inventory status of the proe 
duct they comprise. Examples of this type of product are lubricants, 
paints, glass, plastics, lumber, clothing, gasses, canned foods, dairy 
products, etc, Even though one particular item (or size) may not be 
readily substituted for another, there is still a management need to iknow 
the total quantity and/or value of the basic product for determining pro- 
duction and inventory policies as well as production quantity. For s6xampls, 
if a clothing manufacturer-distributor wanted to produce polka-dot under- 
shorts in various sizes he would need to know the total quantity of cotton 
broadcloth to produce first. (In order to accumulate product data the 
units of an item are converted to the product=-unit by a size factor and 
then sum over all items related to the product. ) 

The simulator can also be used for systems where this feature does 
not occur simply by considering each item as a product and setting ths 
size factor to unity, 

The SRI Inventory Management Simulator/Evaluator is time oriented 
and performs the usual functions of forecasting, ordering, distributing, 
issuing, and record keeping within each time frame of the program, In 
the present form, the simulator has fifty-two time cycles répresénting the 
fifty-two weeks of a year. Within each time frame all the normal inven 
tory functions are carried out in a set sequence which was arbitrarily 


chosen. This arbitrary sequence of functional operations does not bias 





the system since over a long run the system performance will be indepen- 
dent of the actual functional sequence used within a time frame. 

The data for one product and all of its related items are processed 
and massaged completely before beginning simulation of the next product, 
Hence, there is also a cycling cn the various products involved in the 
simulation to produce the overall system data which is desired in the 


summary reports, 


2. Optimal values for system control rules 


The following variables are designed as control variables and form 


a set of parameters which will provide definite policy rules for the 


system: 

RIJK(J ,X) The requisition objective level for the jth item of a 
product at the kth location. The maximzm level in item 
units to which the inventory position is raised at proe 
duction time. 

XIJK(dJ ,K) The reorder level for the jth item of a product at the 
kth location. The levei in item units below which an 
interim branch order is placed, 

Ra) The system-wide reorder level for the jth item of a 
product, The level in item units below which a new pro- 
duction order is placed, 

XI The system wide reorder level for a product, The ievel 
in product units below which a new production order is 
placed, 

XIJK(J,1) The central warehouse reorder level for the jth item of 


a product, The level in item units below which a new 





production order is placed, 


QI The minimum economic production order quantity (E0Q). 


Since the overall objective of the supply system is to minimizes 
its total cost of operation, optimal settings of the control variables 
Could be obtained by determining each of the cost components (i.e., 
holding cost, production ordering cost, and shortage cost) es a *Sunaztioan 
of the control variables and then minimizing with respect to all of the 
variables. However, the dual functioning of the warehouse as a systom 
stock point for interim orders and as a branch outlet unduly complicates 
the mathematics. To avoid this, the warehouse is considered to bs an 
independent branch with demand rate equal to the sum of the direct demand 
rate plus the demand rate due to interim branch orders, This reduces ovr 
two level system to a one level system with independent branches, The 
optimal values for the control variables are then computed according to 
formulas developed in 2, This approximation seems reasonable since 


the interim branch orders are not expected to be too great, 


Economic Production Order Quantity 


The optimal value for the economic order quantity is given in la.) 
po . 
as: QI = { 2/3 1/3 if [ou/t| = 2.25 
1.85 Qw 1 
otherwise 

The "Qw" is the classical or "Wilson" economic order quantity which 
balances the inventory holding costs against the production ordering 
cost. The "f" is a correction factor to account for the demands during 
i 


Wilson, R. H. A New Method of Stock Control. Harvard Business Review, 
Vol. 5, 1926-1927 pp.197=-205. 
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leadtime and the variance of the demand and leadtime distributions, If 
the value of "Q" so determine approaches the shaif life limitations then 
& maximum limit is taken as the demands during a shelf-life period lass 
one and a half times the demands during a production leadtime, This is 
a rather arbitrary setting for the maximum vaiue, No minimum valus is 


set, 


Product Reorder Level 

The optimal reorder level for the product is taken as the differencs 
between the optimal requisition level for the product and the 2conomic 
order quantity. The optimal requisition level for the product is thea 
sum of the optimal levels for the branches, Since reviews of the inven 
tory positions are made at discrete interwais, the time at which the 
reorder levels are reached will fall randomly within the review interval, 
Then the expected demands which occur after the reorder level is actually 
reached and before this fact is recognized will be the demand rate times 
one-half the review period. Likewise, since the average number of units 
demanded per requisition is not necessarily unity the inventory position 
when reviewed will be below the reorder level by one-half of the averages 
number of units per requisition, Hence, the reorder level is to be 
raised by one-half of the demands during a review period plus one-half 
of the expected number of units demanded per requisition to compensats 


for the discrete nature of the review time and demand quantity. 


Branch Requisition Objectives (Other Than Central Warshouse) 
The optimal requisition objective for an item=-location is derived 
as a function of the demand rate, the total product demand rate, the 


variances of the demand and leadtime distributions, the production 
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leadtime, the holding cost and the cost of shortages, The mathematical 
details are given in lad, A maximum and minimum level are sst to guard 
against shelf life limitations and unusual values of the input paraq= 
meters, The maximum requisition objective level is set at the expected 
demands during a shelf life period and the mininum is sét at the reorder 
level plus the demands during a production leadtime and a production 
cycle, 


The formula for the optimal requisition objective level is: 


x 1 Ql ul + 
Rs gic = Pague g gy m1) + Oy Mk Sg By) )/t 


where fijk is demands during leadtime corrected for the variances 
of the demand and leadtime distributions 


Rey is a variance factor 
lL 

1 Q 
(D5 5x -1) is a probability factor 
o? i © of short 

ik s a cost of shortages 
Ks is another variance factor 
Chyk is the holding cost 
Lj is the production leadtime 
t is the production order cycle 


Central Warehouse Requisition Objective 
Like the requisition objective levels for the branches, the optimal 
requisition objective for the central warehouse is determined as a funce 


tion of the demand rate, the variances of the demand and leadtimes 





distributions, the branch reorder leadtime, the holding cost, and the 
cost of shortages, However, in the case of the central warehouse ths 
situation is complicated by the fact that the total demand is the sum 

of the direct demands and the branch interim orders, Since the later 

is a random quantity, it must be computed as an expected value, It is 
assumed that the "tail" of the distribution of the demands at any branch 
can be approximated by a negative exponential whose mean is taken to be 
that of the true distribution, Then the probability of an order occurring 
times the expected size of an interim shipment given that the need occurs 
is the unconditional expected demands on the central warehouse from a 
branch, The details of the computations of the probability of an order 
occurring and the conditional expected size of an interim order are 
given in a. 

The cost of shortages must also be computed as an expected value 
since the cost of shortages pertaining to non-availability of stock to 
meet direct demands is not the same as that for demands arising from 
branch interim orders. For direct demands the cost of shortages is the 
profit lost. The shortage cost applied to the central warehouse for 
not being able to fiil branch interim orders is computed as a function 
of the rate of change of the expected shortage quantity in a branch re= 
order leadtime with respect to the reorder leadtime, As the reorder lead= 
time increases the probability of a branch shortage increases and hence, 
the expected shortage at the branch and resultant costs of shortage will 
increase, The actual cost to be lodged against any warehouse shortage 
is then the expected value of the direct and indirect demand shortage 
costs average with respect to the frequency of occurrence or probability 


distribution of the direct and indirect demands, 
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As with the requisition objective levels for the branches, a 
maximum and minimum levels are established for the warehouse to guard 
against shelf life limitations and unusual settings of the input para= 
meters, The maximum is taken as the total warehouse demands aorine a 
shelf life and the minimum is the sum of the reorder level and the total 
demands during a production order cycle and production order léadtime for 
the warehouse, 

Branch Reorder Level 

Although an optimal reorder level formula is darived in{2.| which 
is analogous to the optimal requisition objective .ormula, it is not 
utilized in the simulator, The branch reorder level is taken to ve the 
demands during a branch interim reorder leadtime, The simulator does 
provide a means of increasing or decreasing this level to study the 


sensitivity of the overall system to changes in the reorder point, 


Central Warehouse and Item Production Reorder Levels 

The central warehouse production reorder level is taken to be ths 
total demands on the warehouse during its production leadtime, 

The system-wide item production reorder level is the demands during 
a branch production leadtime summed over ali branches including the 
central warehouse, 

Provisions are made in the simulator to adjust the reorder level 
settings in order to study the effects of changes in these variables upon 


the overall system performance, 


3. Functional routines of the simulator 
This section will give a general description of each of the funce 


tional routines within the simulator, Figure 1. is a block diagram of 
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these functional routines. 
A. Start New Product 

The time cycle begins at this section and continues through the 
END-OF=WEEK HOUSEKEEPING section. The product, item and branch data 
cards are read into the program at this point. All the parameters of 
the system have now been set and the simulation of a product can now 


commence, 


B. Initial Forecast of Demands 

The twelve monthly demands for each item=-location are averaged to 
provide a mean weekly demand rate, There is some correlation introduced 
here between the predicted demands and the actual demands that the branches 
will experience later since both the initial predicted demands and the 
actual demands are based on the monthly demands, Howsver, it is doubted 
that this correlation is sufficient to appreciably influence the system 


performance except under rather erratic demand patterns, 


C. Initialize Order Queues 
The Order Queues are initialized to ensure that any data left in ths 
queues from the previous product simulation is not erronsously introduced 


into the present product simulation, 


D. Compute Optimal Base Rules 

The computation of the optimal base rules is really the heart of 
the simulation program. The decision rules or policies developed here 
will in effect make or break the system, In the present form of the simue 
lator the optimal requisition objectives and reorder levels are those 


published in bullet The subroutine "RXCOMP" computes these levels based 


LO 





on the initial mean demand rates of the item=-locations, These levels 
are computed first in item units and then converted to weeks of supply, 
In the casé of the central warehouse, its optimal requisition objective 
and reorder level is given as so many weeks of total item demand (all 
branches), Likewise, the system item reorder level is given as so 

many weeks of total item demand. Ths optimal requisition objective for 
all branches other than the central warehouse is given as so many weeks 
of demand at the individual branch, The product requisition level is 
taken as the item=-location requisition objectivés summed over all items 
and all locations, The optimal or economic order quantity (E0Q) is ex- 
pressed as the ratio of the initial economic order quantity to the square 
root of the total product demand rate, (See figure 2, for a block dia- 
gram of this subroutine). 

By expressing the initial requisition objectives, reorder levels 
and economic order quantities in thess ways, they can be corrected or 
updated with the subsequent demand experience, For example, by express] 
ing the requisition objectives for the central warehouse in weeks of 
total item demand, the corresponding number of physical units to use 
for the "current" week's requisition objective will bs the “optimal” 
coverage period times the predicted item demand rate for the current 
week, Hence, as the predicted demand rate changes so will the requisi- 
tion objective, The same principle applies to all the requisition ob-= 
jectives and reorder levels, In updating the economic order quantity 
the square root of the product demand rate is used since the demand rate 


enters the Wilson's E0q" formula as the square root, 


“ Wilson, R. H., A New Method of Stock Control, Harvard Businsss Review, 
Vol. De 1926=1927, Pp. 197=205, 


Li 


E, Predict Demand and Generate Actual Demands 

The subroutine PREDICT is used to forecast the weekly demand rates 
for each item-location, each item, and each product, The forecasting 
scheme is a nine week moving average, There is no trend correction (or 
smoothing) applied to the forecast, The item demand rate is the sum of 
the direct demands (does not include branch interim orders on the central 
warehouse) over all locations and the product rate is the sum of all item 
rates converted to product units, 

Because the input data is expressed in months a complicated indexing 
system is necessary to seiect the most recent "nine week period" on which 
to make the forecast, (See the detailed description of this selection 
process given in the card write up section fer card 137). 

The subroutine GENDEM scales the monthly demand data down to a4 
weekly rate. The indexing used here to select the appropriate month 
for scaling approximates the calendar month-to-week relationship as clos« 
as possible, (See the detailed description of this subroutine given in 
the card write up section for cards 913-939). 

Figure 3, is a block diagram showing the relationship of these sub- 


routines to the main program, 


F,  Gompute R and X Levels 

At the beginning of each time period in the simulation of a product, 
the requisition objectives and reorder levels are updated in accordance 
with the demands predicted for the period. Hence, these policy levels, 
as expressed in physical units, are changing each week but the "optimal" 
period of coverage remains constant. Under the SRI rules these "optimal" 
coverage factors are the outputs of the RXCOMP subroutine while under the 
alternate rules the coverage periods are given as the branch and plant 
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(or central warehouse) goals. It should be noted that under the later 
rules the reorder levels are taken as a percentage of the requisition 
objective, This percentage can be arbitrarily set by changing the value 


of the variable RGX on card 31. 


G. SRI Production Ordering Procedure 

The decision whether or not to place an order for new production 
is made on the basis of the system-wide inventcry position of the product 
and/or of the item, and the inventory position of the item at the Central 
Warehouse, Each product and each item within a product classification 
has an assigned "optimal" system-wide reorder level, In addition, the 
central warehouse has an assigned reorder level for each item in its 
local inventory. Each branch also has an assigned reorder level but the 
fact that a branch (other than the central warehouse) has fallen below 
its item reorder level will not trigger an order for new production, It 
will, however, trigger an order for an interim branch requisition on the 
central warehouse stock, When any one of the thres inventory positions 
(i.e., the system-wide product level, the system-wide item level, or the 
central warehouse item level) falls below its reorder level a decizion 
is made to place an order for new production. The system-wide inventory 
position of the item is reviewed to prevent imbalances between items of 
the same product since as mentioned before no provision is madé to sube 
stitute one item for another. The item inventory position at the central 
warehouse is checked to insure that any imbalance between branches (with 
respect to the same item) does not cause eventual shortages at the branches 
with the deficiency. Since there is no redistribution between branches, 
the central warehouse must maintain its stock levels even though there 
might be sufficient stock of the product and/or the item in the system, 

13 
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Once the decision is made to place a production order it becomes 
necessary to determine what quantity to produce, Since the underlying 
philosophy of the supply system is to anticipates branch needs as accurately 
as possible and to automatically resupply them with as few interim ship- 
ments as possible the production allocations are based on the latest 
possible inventory position, Hence a deferred allocation procedure is 
utilized, The initial production order quantity is taken as the difference 
between the branch requisition objective and inventory position levels 
summed over all branches and all related items of a product in product- 
units, The quantity is then adjusted to the next higher production batch 


size (you might manufacture 10,100 yards of broadcloth, but not 10,001}. 


H, Common Production Ordering Section 

After determining the production order quantity of the product an 
initial (or tentative) allocation is made to the related items and an 
item-quantity allocation is also made to the central warehouse, These 
allocations are made to record the “dues from production" on the central 
warehouse records and the system-wide product and item records, This is 
necessary in order to prevent production orders from being placed apain 
at the next review based on the same needs that trigger the current pro= 
duction order, It should be noted that item-quantity allocations and 
central warehouse allocations are made to the nearest packaged lot (e.g., 
case lots). The production order quantity, delivery date, allocations 
to item-quantities, and allocations to the central warenovse are tinen 
stored in a queue to simulate the production-run leadtime delay, Pro= 
duction order "dues" are not established at the branches but only for the 
central warehouse and the system-wide item and product records, This is 
equivalent to assuming that production order allocations are received a% 
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the branches prior to making the allocations for the next production 


run. 


I. Fill Production Orders 

When the delivery date for a new production order comes due, the 
initial allocations to item-quantities and to the central warehouse are 
canceled and the system item needs are recomputed based on the current 
inventory positions of the branches, At this point in time ths production- 
run component of the production order leadtime has already expired and only 
that fraction due to the waiting and shipping components remain, Since 
the requisition objective was determined as a function of the total lead= 
time, the branches would be over supplied if the deferred allocations 
were based on the optimal requisition quantity. To avoid this situation, 
the requisition objective is adjusted downward by thé amount of the average 
demands during a production=-run leadtime, In determining the needs at the 
central warehouse both its direct demands and branch reordér demands are 
considered. Since the quantity of the basic product produced may not be 
the same as the sum of the current item needs the item ailocations aré 
taken as the current needs times a proportionality factor which is ths 
ratio of the quantity produced to the sum of the current item oeeds, In 
the simulator this quantity is then rounded to the nearest packaged lot 
(case), These quantities are then delivered to the central warehouse 
whose records are updated, Immediately afterwards, allocations are made 
to each item=location based on the adjusted requisition objectives and 
current inventory positions and placed into the branch order queue to 
await delivery at the appropriate time (see figure 5). The allecation 
to the warehouse is in effect the residual quantity not required to fill 


the branch needs, It is tacitly assumed in the simulator that the 
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quantity in stock in the warehouse plus the quantity needs of the 
warehouse are sufficient to cover any incréase in total branch needs dus 
to scaling and rounding. This assumption seems valid except in ths ine 
stance the warehouse does not stock the item and has had no demand and 
even in this case the probability of exceeding the production quantity 


to any significant extend would seem very small. 


Jo Alternate Production Ordering Procedure 

Under this procedure production crders are placed at seach scheduled 
review period. Ths production order quantity is based solely on the 
needs of the central warehouse since under this concept all branch needs 
are filled by placing orders on the central warehouse, No direct allo-= 
cations are made to the branches. Hence, this system can be thought of 
as a "pull" system, <A minimum production order quantity is sst to ensure 


economic plant operations. 


K, Fill Branch Orders 

This section simulates the delivery of both branch interim orders 
and production allocations at the branches, The branch ordsr queue is 
screened to see if the delivery date of any order is due and, if so, ths 
on-hand, on-order, and in-transit records are updated accordingly, 


Figure 4. is a block diagram of this section of the progran. 


Lo Fill Demands on Branches 
In this section the actual weekly demands are presented to the 
branches for issue, If the on-hand balance is nowt sufficient to cover 
the weekly demand, stock may be issued in partial fulfillment of the re- 
quest or the demands may be canceled. If the stock out policy selector 
(ICON3) is set at 1, the requests are partially filled, if it is set av 
Ves 





2, the request is canceled. If the request is only partially satisfied 
the balance is canceled since backorders are not considered in the 

present form of the simulator, The tallies representing the item-location 
on-hand balance, inventory position, satisfied demands, and total demands 
received to date are updated as well as the system, item and product ine 


ventory positions, 


M, SRI Interim Branch Ordering 

Although the system is designed to resupply the branches by allocas 
tions from the production orders in sufficient time to mest their expected 
needs, unusual demands may reduce a branch's item level sufficiently to 
warrant an interim resupply order. Consequently at each review tims the 
item-location inventory positions are screened against their reorder 
points, If the inventory position is at or below this level, a decision 
is made to resupply the branch with the particular item(s) from the 
central warehouse, 

After making the decision to place an interim branch order we must 
determine the quantity to requisition, As in determining the branch nesd 
for the deferred allocation we must adjust the optimal requisition level 
downward by the average demands during the period between the plating of 
the most recently delivered production order and tne present time, In so 
doing we are in effect compensating for only those demands which are aq 
bove the normal rate and hence, are aliowing the system procedures to 
function with the least interference, The difference between this adjusted 
requisition level and the current inventory position is ths “optimal” ine 
terim order quantity, However, in some instances, depending on the rela 
tive lengths of a production order cycle and the production leadtime, 


this quantity may be too low and a minimum quantity must be set, 
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In determining the minimum value of the interim order quantity when 
the production leadtime is short compared to the production order cycle, 
it may be assumed that the need for reordering will occur within a pro= 
duction leadtime, (i.e.) after an order for new production has been 
placed but prior to its receipt at the branch, It may also be assumed 
that this need will occur within the latter half of the production lead-= 
time. Hence, if the branch inventory position is raised to the reorder 
level plus one half of the expected demands during a production leadtime 
this quantity should suffice until the next production allocation is 
received, 

On the other hand, if the production leadtime is long compared to a 
production cycle (period between placing orders for same product) the 
deliveries will be more frequent with perhaps smailer quantities and the 
need to order will occur most likely within a production cycle and pre= 
sumably within the last half of that cycle if the system is performing 
fairly close to expectations, Hence, in this case the minimum interin 
order quantity should be the difference between the reorder level plus 
the expected demands in one half a production cycle less the actual in- 
ventory position. 

The actual interim branch order is then the maximum of the "optimal" 
quantity and the smaller of the minimum order quantities. In the simula] 
tor this order quantity is then taken as the next highest case or pallet 
lot depending upon whether the additional holding costs aré greater than 
the additional handling cost incident to processing broken pallet loads. 

Once the order quantity is finalized, the simulator updates the 
warehouse records and stores the order with its delivery date in the 


Branch Order Queue (see figure 5.). At the appropriate time the order 
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quantity is removed from the queue and the branch records are updated 


accordingly. 


N. Alternate Branch Ordering Procedure 

Since branch orders are only placed at the scheduled review times 
the simulator first checks to see if a review will be made in this time 
frame, When a review is to be conducted, the stock status of the central 
warehouse is first checked to ensure availability of stock, If the cen- 
tral warehouse inventory position is below its reorder point the require] 
ment at each branch is cut in half to supply some stock to all the branches 
which require replenishment, It is implicitiy assumed in both instances 
that the total quantity supplied to the branches does not exceed tne 
warehouse's on=hand balance, It should be noted that there is no provie 
sion made for placing branch orders at times other than at the scheduled 
review date. For example, if a zero balance is reached no action is taken 
until the next review regardless of how long that interval might be, if 
the review interval is one time frame unit there are no difficulties, 
However, if the review period was several time frames in length and if 
performance of the system under the alternate rules were to be compared 
to its performance under the SRI or any other set of rules this limitation 


might be significant. 


O. End-of-Week Housekeeping 

A miscellaneous section where the average inventory on hand, the 
average value of the inventory in-transit and the number of weeks of 
stockout tallies are updated. There is also an optional print statement 
to produce an inventory status report for each item=-location for each 


time period, This section also marks the end of a time frame (week). 
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The program either returns to the predict section to begin the next 
time frame or continues on to print the output reports if the present 


product simulation is completed, 
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SCHEMATIC DIAGRAM OF SIMULATOR 
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| PILL PRODUCTION ORDERS 

Simulates delivery of the item from 
the plant to the warehouse and makes 
final allocations to items and to 
branches, Places allocations in the 
branch order queue. 


| FILL BRANCH ORDERS 


| Simulates the receiving of both 
interim orders and production 
_ allocations at the branches 
} 


"FILL DEMAND ON BRANCHES 


Simulates issues or sales of items 
at the individual branches and up= 
| dates inventory records 


| 

v 
___ SRI INTERIM ERANCH ORDERING 
| 
_ Determines whether a branch needs 
| an interim replenishment and if so 
' the appropriate quantity to order 
from the warehouse 


V/ 


SRI PRODUCTION ORDERING 
Decides whether or not to order and 


___ computes order quantity 
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wee OMMON PRODUCTION ORDERING = 


| Makes initial allocations of the 

| order quantity to the central ware- 
| house and to items for record pur- 
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END OF WEEK HOUSEKEEPING 
Miscellaneous section where certain 
| inventory tallies are updated 
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MAIN PROGRAM 
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SUBROUTINE RXCOMP (ICON) 


Compute economic order quantity, production | 
order cycle time and the requisition rate ! 
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Compute the expected size of branch interim 
orders on the warehouse —s—| 


| Compute the expected cost of shortage at th: 
central warehous 

Compute the system-wide item reorder level 

and the warehouse item reorder level 
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Compute the requisition level for the wars] 
house and the reorder level for the product _ 





PREDICT DEMANDS AND GENERATE ACTUAL DEMANDS 


—————— 


| MAIN PROGRAM | 


Gant Predict | 


Computes the predicted demands for 
the current time frame using a nins 
week moving average 
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Return to 
main {program 









Call GenDem 






SUBROUTINE GENDEM 


Scales the appropriate monthly 
demand to a weekly demand whith 
is used as the sales for the 
present time frame 










Return to 
main |program 






COMPUTE R and X LEVELS 


Updates the optimal requisition 
and reorder levels using the 
current predicted demand rates 








Continues 


Figure 3, 
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RECEIPT OF ORDERS AT BRANCHES 


MAIN PROGRAM 


FILL BRANCH ORDERS Section 


Call Queue B(2) 


S¢treen the 


queue again 
SUBROUTINE QUEUE (N) 







Screens the order queue to determine 
if any orders are due in this tine 
frame 


Orders due 


Updates branch inventory 
records 


GO TO Statement 
/ 
FILL DEMANDS ON BRANCHES 


Continues 


No orders due 


Figure 4 
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PLACE ORDERS IN BRANCH ORDER QUEUE 


MAIN PROGRAM 


FILL PRODUCTION INTERIM BRANCH 
ORDERS ORDERING 








CALL PBO(2) CALL PBO(1) 


SUBROUTINE PBO( ICON) 
Return Updates inventory and cost records Return 


CALL QUEUE B(1) RETURN 


SUBROUTINE QUEUE B(N) 


Stores order until delivery 
date is due 


Figure 5, 
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CHAPTER ITI 


DETAIL DESCRIPTION OF THE FORTRAN LANGUAGE 
CODING OF THE PROGRAM 


Card No. Description 

iL Program Name 

2-4 COMMENT cards 

5-8 Dimension Statement. These variables are used in main 


program in addition to those variables that are listed 
in common, There are 25 possible branches and 20 
possible items or sizes, 

9-27 Common Storage. Dimensioning and placing of variables 
used in main program and subroutines into common 
storage locations, 

28=37 Start of Program. Initialize various variables and 
setting End of File flag to 1, Program stops on End 
of File flag set to 2. 

38-41 Read in the Control, Base Rules parameters, and 
Branch data cards, These decision rules and parameters 
are now set for the entire run. If different decision 
rules and/or parameters are to be used new runs will 
have to be made separately. The variable 7 is used 
as a dummy for both the obsolescent cost and the 
variable VLPl which is no longer used, 

42 The opportunity cost is read-in in the same format as 
the physical storage cost per unit (APB and APP) and 
must be multiplied by .1 to obtain correct valus, 

This is done solely for convenience of formatting. 
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43-47 


48-50 


als 


De. 


D3-94 


59 


56-68 
69-72 


Control Data, Decision Rules parameters and Branch 
data are printed out to check correctness of read-in 
and also to have this information listed on ths cute 
put sheets for future identification and reference, 
Leadtime components are readin in "days" and ars 
converted here to "weeks", 

Call SUBROUTINE RXCOMP to initializes variabiese used 
for system totals, 

Call SUBROUTINE RDTP to read=in the first produst 
characteristics, 1.e., the base number assigned to 

the products the shelf life expressed in months; the 
branch goal expressed in days; the plant goal (central 
whse, stock) expressed in days; the production lead- 
time expressed in days; the batch size increment ex- 
pressed in product units; minimum batch size «expressed 
in product units. All other product cards are read 

by card number 105. 

Use a do=-loop to set all of the "intransit" leadtimes 
for the branches to same value. This is saquivalent 

to assuming branches are equidistant for shipping 
purposes from the main warehouse/plant. 

Re=-set the warehouse shipping leadtime to zerc. Since 
the warehouses is located at the manufacturing plant 
there is no shipping time involved, 

Use a do=loop to initialize branch variables, 

Start New Product: Set the time clock to zero and 


initislizs the product variables, 
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73-79 


80 


81-82 


83 
84-95 


Convert shelf life, branch goal, plant goal, and 
production leadtime to "weeks", and transfer the 
values cf batch size increment (BSI), batch size 
mean (BSM), and base number (BASE) to permanent 
storage locations, It should be noted that all 
product, item, and branch characteristic data ar« 
first read into the AO array and then transferred 

to permanent variable name locations, When ths 
program was used by SRI this data was read-in from 

a tape rather than cards and having a standard format 
for presenting information simplified the problem of 
handling massive amounts of data, 

Call the subroutine RDTF for the first item of the 
present product, Here the size code, item number, 
cartons per pallet, units per cartons, weight per 
unit, production unit cost, selling price, and size 
conversion factor (# of product units per item unit) 
are read. All subsequent item, and branch data cards 
related to the present product are read by the "Call 
RDTP" on card number 105. 

Set the count for the different items (sizes), As 
each item card is read the program returns to card 
82 where the item count is increased by one. 
Initialize the inventory position of the item. 
Initialize the branches! (item-location) inventory 
records; satisfied demands, order count, stock out 


weeks, stock out counts, broken case count, broken 
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96-104 


105=124 


Wes) 


pallet count, on order quantity, average inventory 
position, total demands and the branch flag. 
The product number (four digits) and the size code 
(two digits) are listed as a six digit number on the 
data card, Cards 96, 97, and 98 ars used to separate 
these two numbers and store them in permanent loca- 
tion and to scale data as necessary, 
On the first pass through this section of the pro- 
gram we have already read the first product card and 
the first subsidary item card, When we call "RDTP" 
the first time at card 105 we will read a branch data 
card which will leave a value of 3 in the AO(1) loca] 
tion, By testing this variable, the next two cards 
will send us to card 110 where we select out the 
proper value of k by matching the branch code array 
to the branch code number read from the data card, 
Having now determined the proper item and branch 
numbers (j and k) we transfer the input data to perma- 
nent storage arrays. As a bonus operation we also 
compute the mean item=-location demand rate and the 
item and product inventory positions, 
This card sents us back to card 105 where we call 
RDTP again and read the next card. Four sequences are 
possible at this juncture: 

(1) If the next card is another branch usage 

data card we will repeat the same sequence 


as outlined above, 


31 





(2) 


(3) 


(4) 


If the next card read is an item record 
ecard, AO(1) will be set to 2,0 and this 
time card 107 sends us back up to card 8&2, 
In this way the machine recognizes that no 
more branches are carrying the previous 
item and that a new item is now being cone 
sidered (but still under the same product). 
The item count, NJ, is increased by 1 and 
we initialize all the item and branch 
variables for the second item, We are then 
back to card 105 where we call RDTP in 
search of branch usage data cards for the 
new item (size), We now follow sequence 
(1) again, 

If a product card is read, AO\1) will be 
set to 1.0 and card 107 will send us to 
card 129 where we commence simulating for 
the previous product and subsidary items, 
Hence, the data for the next product ré- 
mains in the AO array until the simulation 
is completed for the previous product, at 
which time the program returns to card 69, 
If, when card 125 sends us to card 105, a 
blank data card is read, the RDTP subroutine 
sets AO(1) to -1,.0 before returning to the 
main program, Hence, when we test AQ(1) 


this time we are sent to card 108, which 
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126-130 


131 


132 


sets End of File to 2, We are then sent 

to card 129, We have completed reading in 

all the data cards and will simulate for 

the last product, print out the summary 

statements and stop the program, 
Initialize the Order Queues. Before beginning the 
simulation we wish to initialize the Production Order 
Queve and the Branch Order Queus. The present routine 
sets the initial values to 9.0 so in effect we ars 
simulating a system from its beginning, If this is 
not desirable the program could be modified so that 
representative values could be read in to simulate a 
system in progress. 
Compute R and X Levels, Hers we check to see what 
set of Decision rules will be used in controlling 
the production ordering and branch shipments, If 
ICON] has a value of 1 we use the SRI rules, other= 
wise we use the alternate set of policy rules. 
The "RXCOMP" subroutine is called to obtain the 
optimal policy rules for the present product, The 
requisition objectives and reorder levels for seach 
item location, the system-wide item and product 
reorder levels and the minimum economic order quantity 
are computed within this subroutine, The requisition 
and reorder levels are expressed in weeks ot supply 
rather than in item or product units so that they 


may be varied as the current demands varies, 
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aia 


= 135 "Generate Actual Demand, Predict Demand, and Compute 
X and R Levels" heading, 

136 Start ciock by setting time equal to 1. 

ai? Call SUBROUTINE PREDICT to make forecast of the 
weekly demand rate for each item-location, each 
item, and for the product. Demands are computed 
by taking a weighted average of the demands in 
those months in which the previous nine weeks 
occurred, For example, the prediction for the 
first week is based on the demands occurring in 
the nine weeks (the previous two months) immediately 
preceding the simulation period, It is assumed that 
five weeks occurred in one of the months and four 
weeks in the other, Hence, the demands in the past 
two months are weighted by five=ninthe and four= 
ninths, This weighted monthly average is then 
scaled by 7/30 to obtain a weekly average, In the 
Simulator, the demands of the eleventh and twelveth 
months are also used as the demands for the two months 
prior to the simulation period, 
In order to select the appropriate months in which 
the past nine weeks occurred, a complicatsd indexing 
system is needed. The schedule below indicates the 
results of this system for the first three time 
frames, The Jl index provides 4 means of assigning 
the nine week period to the appropriate months as 


the time index moves through the 52 simulation periods. 
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138 


139 


140 


The DM(i) is the demands in the i-th month. 


Jl Index Monthly Demand 


45 DM(11) 
46 DM(11) > 
47 DM(11) | 
48 DM(11) | 
Bi 
49 DM(12) 
50 DM(12) B 
51 DM(12) T 
52 DM(12) : 
53 DM(12) 
54 DM(1) 
55 DM(1) 


e@ 5S. 
Call SUBROUTINE GENDEM to generate the actual weekly 
demands by scaling the monthly demands inputtsd from 
the data cards by 7/30 and rounding to the nearest 
whole unit. Note that the predicted demands are de= 
signed as DIJK whereas actual demands are designated 
as DWIJK. 
Branch to SRI rules or the alternate set as deter= 
mined by the value of the inventory rule selector 
(ICON1) ¢ 
Compute the economic order quantity (E0Q) for the 
current week, The variable FLQI is the ratio of the 
optimal economic order quantity to the square root of 


the mean product demand rate, By multiplying this 


By, 





141 


142-143 


144-145 


146 


147 


148-150 


ratio by the current predicted product demand rate 

we are varying the EQQ to reflect the present demand 
trends. The variable FLDI is the maximum EOQ due to 
shelf life limitations, 

Compute the production order cycle as the E0Q coverage 
period plus one half of the review interval. 
Initialize the product requisition and set up a 
do-loop on items, 

Compute the requisition levels for the central ware= 
house. The variable XRIJZ represents the optimal 
coverage period (of total item demands) and is con- 
verted to physical units by multiplying by the cur- 
rent predicted demands for the item. The variable 
XDIJK represents the maximum level due to shelf life 
limitations. The minimum level is given as the re= 
order point plus the demands in a production order 
cycle and leadtime. 

The reorder level for the central warehouse is taken 
as the total demands (including interim orders) during 
@ production leadtime. The variable XXIJZ is the ratio 
of the mean total demand during a leadtime to the mean 
direct demand rate. 

The system-wide item reorder level is also taken as 
the demands during a leadtime, 

Initialize the item requisition variable, establish 

a do-loop on locations and screen the branches to ses 


if the item is carried, 
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Hee 153 


154 


DD 
155A-156 


J 


158-159 


160 


161-167 


The item-location reorder level is taken as the 
demands in a branch interim order leadtime, 

The optimal requisition objective for the individual 
branches are computed by a formula analogous to that 
used for the central warehouse. However, the pre- 
dicted weekly demands at the branch and its produc- 
tion order leadtime are used rather than the total 
item demand rate. (See cards 144 and 145) 

Sum the requisitions objectives for each item over all 
branches except the central warehouse, 

Continue statement for do-loop on locations, 

Compute the total requisition level for each item 

and each product, 

The product reorder level is the product requisition 
level less the economic order quantity. 

The variable FLPDIM represents a specific time, It 

is the time at which the most recently delivered 
production order was placed, The difference between 
this time and the present time represents that portion 
of the resupply cycle which has already expired. This 
period is used to adjust the branch requisition level 
when computing the order quantity for interim ship= 
ments, (See cards 250, 251) 

Jump statement to skip over the following section 
which is used for the alternate policy rules, 

The requisition objectives and reorder points ars 


calculated for the central warehouse and branches 


af 


168-170 


171-172 


iS 


174-179 


180 


under the alternate policy rules. The warehouse 
levels are a function of the total item demand and 

@ predetermined goal, The branch levels are func- 
tions of the local demand rate for a given item, 

FILL PRODUCTION ORDERS Heading 

Screen the production order queue to determine if 

any orders have been completed by the plant. If 

there is a positive order value due at this time, 

the program will continue on to simulate receiving 

the order in the central warehouse and making 
allocations to the individual branches, 

Check to see which set of policy rules are in use. 
Under the SRI rules the program will compute the 
current needs of the branches and make automatic 
shipments to satisfy these needs. Under the alternate 
rules all production orders are delivered to the ware= 
house for storage. The branches place their orders 

on the warehouse as their reorder points are reached. 
When production orders are initiated a temporary 
allocation is made for each item of the product basad 
on the system needs at that time, and the system item 
and product records as well as the central warehouse's 
records are adjusted for the production "dues", These 
initial "dues" are cancelled here prior to making the 
actual allocations to items, 

Compute the warehouse needs for each item of the 


product. 
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181-184 


185-186 


187-189 


190-196 


M7 


198=204 


205 


206=207 


208-210 


Compute the branch needs for each item and sum over 
all branch needs to give the totai allocation of 

the product to each iten. 

Compute the total system current needs in product units 
and ensure that it is a positive quantity, 

The item needs are scaled to the actual quantity 
produced (POBS) and rounded to the nearest case lot. 
The appropriate inventory records are updated to re= 
flect the delivery of the new items into the central 
warehouse, 

Check to see which set of control rules are in use. 
Compute the order quantity and delivery date for 
each item-location other than the central warehouse, 
The order quantity is rounded to a full case lot, 
The value of inventory in transit is also updated 
here. 

The PBO (Place Branch Orders) subroutine is called to 
update the inventory records and cost tallies for 
each item-location, The parameter value of 2 in the 
Call indicates that this is a production allocation 
and not an interim branch order. The PBO subroutins 
then calls the Queue B subroutine and places the 
orders in the queue to simulate the shipping delay, 
After completing the first production order the pro- 
gram returns to card 17] to screen the queue for any 
other orders that might be due. 


FILL BRANCH ORDERS Heading 
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211-215 


216 


217-219 


220-221 


222 


223-224 


225-241 


Screen the branch order queue to determine if any 
orders are due for receipt at the branch, If an 

order delivery date is due the on-hand, on order 

and in-transit inventory records are adjusted 
accordingly. 

This GO TO statement cycles the program back to the 
CALL QUEUE B card to have another look at the queue, 
FILL DEMANDS ON BRANCHES heeding 

Set up do=-loops on items and locations, 

Check to see if Branch carries this item, If not, 
machine jumps to the next K value and checks again, 
Check to see if there has been any demands this week 
for the item-location., If so, the total demands to 
date tally is increased by the current demands, 

If the on hand balance is positive we compute ths new 
on hand balance (TP) and the quantity of demands 
satisfied. If the new on hand is also positive we 
jump to statement 523 to update the branch, item and 
product inventory records, If the new on hand balance 
TP had been negative or zero the program checks to se 
which issue policy is in use, If issues are to be 
made to a zero balance (a setting of 1) the old on 
hand balance becomes the satisfied demand (TP2) and 
the new on hand balance is zero, In this case the 
stockout count tally is also increased by 1 since only 
partial issue could be made, If on the other hand the 


policy is not to make partial issues the old on hand 
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242 
243-245 
246-248 


249 


250 


251-252 


balance is placed in TP and the satisfied demand (TP2) 
variable is set to zero and stock records are updated 
accordingly. 

Branch to the SRI rules or alternate policy rules, 

SRI BRANCH ORDERING PROCEDURE heading 

Set up do=loops on items and locations and test to 
see which branches carries the item, 

Review the inventory position to determine if a branch 
interim order should be placed, If the inventory 
position is at or below the reorder level an order 
will be placed, 

Compute the time period between the present time and 
the time at which the most recently received produce- 
tion order was placed, 

When placing a branch interim reorder the branch ine 
ventory position is brought up to its optimal requisi-= 
tion level, To do so would overstock the branch and 
interfere appreciably with the automatic production 
allocation system, The policy is to supply the 

branch with enough stock so that at the time just 
prior to receipt of the next production allocation 

the branch will be at its reorder leve), Hence, the 
optimal requisition level is adjusted downward by the 
expected demands during the period discussed on card 
249. A minimum reorder quantity is taken as the re-= 
order level plus the demands in half a production 


leadtime (or cycle whichever is less) minus the 


253=254 


ea) 


256-257 


258-259 


260-261 


262 


263 
264 


current inventory position. 

Check to see if rules call for rounding order 
quantity to the nearest carton and round if they do, 
Check to see if the rules call for rounding to pallets, 
Compute cost of issuing broken pallet load for this 
order, The cost is the number of units in the broken 
load times the cost of handling such units (cBPC), 
Compute the extra holding costs that will be incurred 
at branch if a full pallet load is issued. 

If the holding cost is greater or equal to the broken 
pallet handling cost then the order quantity is 

taken in full case lots, If not, the order quantity 
is taken as the next higher full pallet quantity. 
Compute delivery date for the order, Delivery date 
is specified as the time that the last scheduled de— 
livery was made from central warehouse plus the 
scheduled waiting period plus the in-transit time, 
Update the value of the inventory in-transit, 

Call subroutine Place Branch Order to simulate the 
updating of the branch and central warehouse inventory 
positions, the on order quantity, the records for the 
branch, the on hand balance for the central warehouse, 
and the branch cost records, Subroutine Quéue B is 
called from within the PBO subroutine to place the 
order and delivery date in the Queue until the 
appropriate time to simulate receiving the order at 


the branch, 





265 CONTINUE Statement 

266-268 SRI PRODUCTION ORDERING PROCEDURE heading 

269-284 Check Product Inventory position against product 
trigger level, If inventory is above trigger level 
then central warehouse inventory is checked against 
its trigger level. If it is above the trigger level, 
then the item inventory position is checked against 
its reorder level, If all three are above their 
trigger levels no production order is placed, Howe 
ever, if anyone of the three is below its trigger 
level the production order quantity is computed by 
summing the difference between the requisition ob= 
jective and the present inventory position over all 
branches and all items, rounded to the nearest batch 
size increment, The individual item needs for the 
system POAS(J) and the central warehouse POAW(J) are 
retained in order to make the necessary "dues" on 
the item and central warchouse inventory records, 
This eliminates the placirg of another production 
order later on to fill the same needs, 

286=287 The production order quantity is rounded to the next 
higher batch increment. The program now jumps to 
card 319 to complete the ordering routine, The inter 
mediate sections pertain to the alternate policy rules. 

288=318 (See Alternate Branch Ordering and Production Ordering 
Procedures at ths end of this section). 


319=327 Compute the scaling factor due to rounding to nearest 





328-329 


330-332 
333-339 


340=345 


346 
347 


batch size increment (BSI) and the order delivery 
date, Production order allocations to size and to 
the warehouse are scaled to the production quantity: 
the warehouse's on order and inventory position, as 
well as the product and item inventory positions, 

are updated, 

A tally is kept on the number of production runs made 
during the simulation, Finally, the order quantity 
and delivery date are placed in the queus to simulate 
the production delay. 

END-OF=WEEK HOUSEKEEPING heading 

Set up do=loope on items and locations and select 

out those locations which stock the item. Th average 
inventory on hand and stock out weeke records are up- 
dated, 

This is an optional PRIN” statement used mainly for 
checking out the simulator, It provides 4 weskly 
status report for each item=location showing the on 
hand balance, the on order quantity and the inventery 
position. It also provides the system inventory 
position for the produc% and each related iten, 
Average Value of inventory in transit is computed, 
Check the time to see if the simulation for this 
product is complete, If it is, the program will 
print out the item and product summaries, If not, 

jt will return to the Generate Demand Section, set 


the time to the next increment and go through 


Ah, 





348-350 
Bo 
352-358 


g59 
360-386 


387-429 


program again. 

Print Routine for Item Summary, 

Update the page number tally. 

Set up do=loops and select out Branches carrying 
each item and print item=-location data: 

DWIJK Actual weekly demands in units 

RIJK Requisition level in units 

XIJK Trigger (reorder) level in units 

XRIJK Requisition level in weeks 

XDIJK Demands in a shelf life, item-locaticn 
XDIJ Ratio of total warehouse demands for an 


item to total system direct demands for 


that item 
XXIJ System item reorder level in weeks 
XRIJZ Warehouse item requisition level in weeks 


of total system direct demands for that item 
XXRIJZ Warehouse item reorder level in weeks of 
total system direct demands for that item 
XK2 1K Interim branch order leadtime 
SIJK Inventory position of item-location 
Print page heading 
Initialize product and item print variables prior to 
their computation by summarizing item-lcocation data, 
In this section the system performance data is sum-= 
marized for each item with individual branch sub= 
totals. The inventory values are given in both units 


and dollar figures, The effectiveness measurés such 
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430-445 


446-450 


451-463 


464 


465-467 
468-473 


474-486 


487-510 


as the stock turn ratio and the ratio of lost demand 
to total demand are based on units rather than dollars, 
The holding cost, stockout costs, and costs of handling 
broken carton and pallet lots are also computed for 
the summary report, The print statement causes the 
item-location data to be printed, 

The inventory data is summarized by product totals in 
this section, The print statement here causes the 
item totals to be printed on the summary report, 

The product effectiveness measures (stock turn ratio, 
loss demands to total demands, and weeks per stockout) 
are computed, The print statement prints the product 
totals on the summary report. 

The overall system totals (all products) are accumu=- 
lated here for the Item=Product Summary Report, 

Test to see if the simulation of all products is 
complete, If not, the program returns to the START 
NEW PRODUCT section and commences simulation of the 
next product, 

PRINT BRANCH SUMMARY AND STOP heading 

The overall system effectiveness measures are computed 
and the system totals are printed for the summary re- 
port. 

Initialize the system variables prior to recomputing 
by summarizing on items and products for the branch 


totals. 


Compute the system totals by summarizing the branch 
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582 
583-588 
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607-610 


611-612 


totals. The individual branch totals are also 
printed within this do~-loop. 

Print the system totals. 

The cost of all production set ups. 

The total branch variable costs for the system, 

Page number tally. 

Print the page heading. 

Branch on the control rules. If the SRI rules are 
being used, the program will print a Base Rules Sum-= 
mary which compares the system performance under the 
SRI rules with those of the alternate policy set, 
Print branch variable cost factors and the value of 
the average inventory in transit. 

Stop the program card, 


Formats for READ and PRINT statements, 
SUBROUTINE RXCOMP (ICON) 


Subroutine defined with call parameter ICON. 
Dimension the local variables. 

Dimension and place in Common those variables used in 
one or more subroutines of the program, 

Function VARDF (D,C) defined. See SRI TM -3 for 
information on equation used, It computes the rel= 
variance (ratio of the variance to the square of the 
mean) of the demand distribution, 

Computed GO TO Statement for branching into one of the 


three subdivisions of the program, The first is merely 
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637 
638 


639 


640 


641 


an initializing routine; the second is the computaticn 
of the base rules themselves; and the third is the 
recap and printing of management data taken fron the 
base rules. 

Routine for initializing system management variables, 
COMPUTE BASE RULES heading 

Initialize variables to be computed, 

Computation of expected leadtime between reaching a 
trigger level and availability of item for distribu- 
tion or the production leadtime for central warehouse, 
Sum of half of the review interval, the communication 
leadtime, and the production-run leadtime, 

Set L. ©. Up do-loop for locations 

Compute the mean production order leadtime (reorder- 
to-receipt) for each branch. 

Compute the mean branch interim order leadtime 
(reorder-to-receipt) for each branch. It does not 
include the communication leadtime or the production-= 
run leadtime. 

The item-location reorder point is taken as the demands 
during a branch interim order leadtims, Here the 
branch interim reorder leadtime (FLIK) is multiplied 
by the variable FK2, For manegement purposes, studies 
can be made of how the reorder level influences the 
overall performance of the system by changing tha 
value of this leadtime through the multiplier FK2, 


T5 is an intermediate computation in the variance of 
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643 


644 


645 
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leadtime, It is assumed that the times when produc= 
tion is actually needed are distributed uniformly 
between times of review, and that the times of availe- 
ability are distributed uniformly between times of 
regular (scheduled) branch shipments. If the above 
assumption is true, the variance of the reaction com 
ponent of leadtime is PLRI* /12 and the variance of 
"waiting for shipment" time is FSK(K)* a2. 

The rel-variance (ratic of the variance to the sauare 
of the mean) of the production order leadtime is cone 
puted as the sum of the variances of the five indepen-= 
dent components (reaction time, communication time, 
production time, waiting time, and shipping time) 
divided by the square of the mean leadtime, FLIK(K). 
The rel-variance of the production-run component is 

an input parameter and is used here to compute the 
variance of the production-run leadtime., The shipping 
times and the communication times are taken as nearly 
constant with zero variance, 

Tl - Intermediate computational factor used to simplify 
calculations later, 

Set up do=loop on items. This loop is inside the do= 
loop on locations, 

Test to see if branch carries item, 

Dealer or Branch unit cost of an item, Sum the pro= 
duction unit cost CJ(J), unit freight cost, and unit 


receiving cost. 
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649-051 


652 


653 


654 


655 


Compute the rel-variance of the demand rate as a 
function of weekly demand rate and dealer's unit 

cost. (See [1.| pages 21 and 22 for reference) 

DIPR - Rate of arrivais of orders for the product 

(per unit time) computed as the reciprocal of the 
rel-variance of demand. (See (2. page 25) 

FIJK represents the demands during a production 

order leadtime adjusted for the variances of the 
demand and leadtime distributions, (See 1. page 
20), T2 is an intermediate factor in computing the 
variance multiplier (FLKIJK). 

F is the sum of all the FIJK'S converted to product 
units. It is a factor in the computation of the 
economic order quantity for the product, 

Unit holding cost per unit of time is taken as the 
sum of the opportunity cost and physical holding 

cost. (Cost of holding one unit of an item for one 
year at a given branch). 

DICHI - Factor in computing the Wilson economic order 
quantity (E0Q). In order to determine the average 
unit holding cost (over all branches and items) a 
weighted average of branch holding cost is taken, The 
weighting factor is the ratio of the item-location de- 
mand rate to total product demand rate. The quantity 
DICHI is this weighted average times the total product 
demands rate. 


Sum the predicted weekly demands by location to give 
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658-660 


661-662 
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the predicted weekly demand rate for an item, 

Sum the predicted weekly demand rate by items to 

give the product demand rate in the product=equiva= 
lent units. SJ(J) is the ratio of product units to 
item unit. 

CONTINUE statement 

This PRINT statement is used for checkout purpeses 
only. The values of the print variabies are checked 
against those computed manually to verify the accuracy 
of the program, These cards may be removed after the 
program is checked out. 

Computation of Economic Order Quantity heading, The 
optimal economic order quantity is computed as a 


modified "Wilson" EQQ using the following formula: 





Qvtf 2 
Qw 
Q= for |7# = 262d 
21/3 
1.85(Qw f) 
otherwise 
where 
Qv = V2 AD/C,, 
A is set up cost 
D is product demand rate 
Gy, is average holding cost per unit per year 
f is demands during a branch production 


leadtime adjusted for variances, 


Compute the square of f. 
Compute the average holding cost over all branches. 


Compute the square of the Wilson's E0Q divided by 
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669-670 


671-673 


674 


675 


f squared, The variables CM2 and FMU are multipliers 
and would normally be set to 1.0 and 0.0 respectively. 
Test the value of SQI to determine which formula is 
appropriate. 

If we are in the first region of the formula, we 
convert the test quantity, (Qw/t)*, to the appropriate 
formula by taking the square root, adding one and 
multiplying through by f, 

If the second formula is to be used then the square 

of the Qw/f ratio as computed on card 665 is multiplied 
by the constant 1.85 and raised to the one third power, 
This value is then multiplied by the factor f to ob-= 
tain the correct formulation, 

The maximum or upper limit on the EOQ is taken to be 
the demands during the shelf-life less 1,5 times the 
expected demands during the central warehouse's lead-= 
time. The minimum of the optimal and maximum values 
is taken as the final EOQ for the product. This 
quantity is rounded to the next higher unit. 

The discrete adjustment period is the expected re=- 
action time (FLRI/2) plus the expected time betweer 
orders (one half of the reciprocal of the requisition 
rate). 

The adjusted EOQ coverage time is taken as the pro=- 
duction order cycle (the EOQQ divided by the demand 
rate) plus the discrete adjustment period, This 


later period accounts for the fact that demands arrive 
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677-678 
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680 


681-685 
686-693 


694 
695 


696 


in "lumps" and also for the delay between reaching 
a reorder level and taking action to reorder, 

QIPR is the number of orders corresponding to the 
adjusted EQQ quantity. 

PRINT statement for checkout purposes only. 

FLQI is the ratio of the optimal EQQ +o the square 
root of the predicted demands for the product at 

the beginning of the simulation period, This ratio 
is multiplied by the square root of the current pre= 
dicted demand at the beginning of each week in the 
main program to update the EOQQ value, 

FK4 and FK1 are the same multiplier. The first is 
used in the RXCOMP subroutine and the later is used 
in the main program. They are both used to vary the 
leadtime demands when computing the minimum requisition 
level. (See card 720 and 153) 

Initialize the product variables. 

Set up a do-loop on items and initialize the item 
variables which are to be computed by summing cover 
locations within the do-loop on branches. Note thaw 
the central warehouse is included in the sum. 

Test to see if the branch carries the item, 

Compute the dealer's (branch) unit cost on an item 
which included the freight and receiving cost per 
unit in addition to the production unit cost. 
Compute the rel-variance of the demand rate on an 


item location as a function of the predicted demand 


oe, 








697-700 


POL 


702 


TOL 
705-706 


TOT 


708 


203-/il 


T12 


rate and the branch unit co8. 

Compute the demands during a productson lead time 
adjusted for the variances due to leadtine and de 
mand distributions, 

Computs the hoiding cost for #azh ctem at sath Joey 
tion. 

FKIJK is a variance fasioy uséd in competing ube 
requisitioning level for eatn item location, 

The cost of & shortage at a branch ie taken as the 
profit lost when cut of stoek, 

The probability of ordering, 

Compute the optimal requisition objective level fer 
each item-location, (Ses Chapter II for details of 
the formula) 

Any system and/or branch constrain®2 are applied 
here. In utilizing thes constraint factors the 
immediate effect is tc rsduce the requisition ob- 
jective, This, of course will result in a lower 
average inventory and lowar holding costs. Ths 
system constraint factor may be used with or withont 
any branch constraints. 

Compare the optimal requisition level to demands 
during a shelf-life and take the mininum, 

Print the item-losation dscision rules and input 
data for check out purposes only. 

T5 is an intermediate factor in computing the rel= 


variance of leadtime, (See card 641} 
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718 
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720-721 


722 
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726 


Compute the rel-variance of the branch interim 

order leadtime,. 

Compute the demands during a branch interim order 
leadtime adjusted for the variances of the demand 
and leadtime distributions. It should be noticed 
that the factors computed here on cards 712=717 ars 
analogous to those computed for the production order 
leadtime on cards 641-642 and 697-700, These vari- 
ables are used in computing the warehouse's optimal 
Pequisition PeyeL.. 

CKZJK is an intermediate factor used in computing 
warehouse requisition level, 

The branch reorder level is taken as the interim 
order leadtime demands. 

The optimal requisition level is compared to a 
minimum value which is taken as the reorder point 
plus the demands in a production order cycle and 
leadtime. 

Sum the item-location requisition levels over all 
branches other than the central warehouse, 

DZIJK is the mean of a negative exponential distribu- 
tion which is used to approximate the "tail" of the 
actual demand distribution. (See betel pages 29-31) 
FRIJK is the probability that an interim order for 
the (J,K)th item=location will occur following the 
arrival of a production allocation and before ithe 


arrival of the next allocation. Or, the 
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736-737 


p [D in (t+#PLIK) = RIJK ~ yIIK] . 

Compute the expected size of an interim branch 

order given that an order occurs, 

DPR(J) is the expected size (number of units) of an 
interim order for the j=th item, This quantity 
represents that portion of the total demands on the 
central warehouse by interim orders from the branches. 
Compute the inter arrival time of branch interin 
orders, i.s., the number of weeks per interiz order, 
Compute the probability of an order for the j-th 
item. 

Compute the optimal requisition objective level ex-= 
pressed in weeks, (This card would morse logically 
follow after card 722). In the main program, a new 
optimal requisition objective level is computed each 
weeks using this optimal coverage time, The new 
requisition objective level in item units is this 
coverage factor times the current predicted demands, 
The weekly predicted demands are summed over all 
branches except the central warehouse, 

The demands in a production leadtime are summed 

over all branches except the central warehouse, 

T2 is an intermediate computation to be used in 
computing the average cost of a central warehouse 
shortage. 

CSZ(J) is the average cost of a central warchouse 


shortage of the j-th item. 


56 





738 


739=740 


TAL 
742 


743 


ThA 


745 


746 
Th7 


748 


749 


750 


751 


XDIJK(J,K) is the demand in a shelf-life period for 
the j-Kth item-location, 

PRINT statement for checkout purposes only. 

CONTINUE statement for the do=loop on locations, 
Compute the total warehouse demand rate as the sum of 
the direct demand rate and the branch interim order 
demand rate, Prior to this statement DPR(J) was ziven 
as the number of units demanded in a production order 
cycle, Here this quantity is divided by the weeks per 
cycle (STI) to give the interim order demand rate, 
Initialize the optimal requisition objective (expressed 
as weeks of supply). 

XDIJK is the demand of an item over a shelf lifs 
period at the central warehouse. This quantity is 
used as the upper limit on the requisition objective, 
XK2IK(1) is the reorder time for the central ware-= 
house times a multiplier. 

Compute the total direct demands per week for an itsm, 
Compute the total direct demands in a production lead= 
time for an item. 

Compute the total direct demands in a production lead= 
time for a product in product equivalent units. 
Compute the total direct demands in a production lead~ 
time for a product considering branches only. 

Compute the total direct demands per week for a pro= 
duct considering branches only. 


The system item-reorder level is computed as system 
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756 
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758 


759 


demands in a leadtime. The multiplier FK5 can be 
varied to increase or decrease this quantity. 

Compute the total direct demands for all items of 

One product placed on the central warehouse, 

XXIJ(J) is the system item reorder level expressed 

in weeks, In the main program this coverage period 
times the current predicted demand rate is used as 

the optimal reorder level for the current time period, 
XDIJ(J) is the optimal reorder level for the central 
warehouse expressed in weeks of total demands for the 
item in the system. It is used in the main program 
to compute the optimal reorder level in physical units 
by multiplying it by the current item demand rate. 
XJ1(J) is the central warehouse reorder level for the 
j-th item taken as demands in a production leadtime. 
XXIJ2Z(J) is the central warehouse’s reorder level 
expressed in weeks of total item demands in the 
system, 

Sum the item requisition objectives to determine the 
product requisition objective. 

CIJl is the unit cost of an item at the central ware= 
house; production cost plus delivery and receiving 
costs. The later costs are given in dollars per hun- 
dred weight and must be scaled down by .Ol1 to give 
dollars per pound, The weight factor, WJ(J), is given 
in pounds per unit. 


Transfer the unit cost figures to their permanent 
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storage locations. 

Compute the rel-variance of the direct demands on 

the central warehouse, and transfer this quantity to 
a permanent storage location, 

CHJJ(J) is the holding cost per unit per week, AZR 
is the opportunity or interest cost expressed in 
dollars per dollar of inventory per year and APP is 
the physical holding cost expressed in dollars per 
product unit per year. 

VDPR(J) is the rel-variance of the total demand oa 
the central warehouse for each item. 

DIPRPR is the requisition rate corresponding to the 
total demands on the central warehouse for each item 
summed over all items to give requisition rate for 
the product, 

PRINT statement for checkout purposes only. 

CONTINUE statement for do=loop on items, 

FKA is a variable multiplier used to increase or de= 
crease the leadtime-demands component of the mininum 
requisition objective. It is used in the RXCOMP sub- 
routine in determining “optimai" values of the requisi- 
tion objective of item for all branches including the 
central warehouse, In the main program, however, FKL 
is used when updating the item requisition objectives 
for the branches and FK3 is used when updating the 
item requisition objectives for the central warehouse, 


This provides a means of using a different multiplier 
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on the central warehouse than on the other branches 

if so desired. 

Sum the product, requisition rate at the central ware= 
house to give a system requisition rate at the central 
warehouse, 

Transfer the total of branch requisition objectives 
(except central warehouse's) to another storage loca- 
tion for later summarization for system total. 

Compute the number of branch interim reorders per 
production order cycle. 

Compute the values of little f and k, Capital K is 
used in determining the optimal requisition objective, 
Compute the expected value of the cost of a shortage 

at the central warehouse, Since shortage at the 
central warehouse may mean a failure to meet direct 
demand and/or a branch interim order the average value 
must be taken, This cost factor is further complicated 
since a failure to meet an interim order may or may not 
result in an actual shortage at the branch prior to the 
receipt of the next production allocation, 

BPIJK ig a computational factor representing the prob- 
ability of ordering. 

RPJK is the value of the optimal requisition objective, 
Modify the optimal requisition objective value for 
system and warehouse constraints, 

Compare the optimal value of requisition objective to 


total warehouse demands in a sheif-life period and take 
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minimum quantity. 

Compare the optimal value of requisition objective 

to the lower limit and take the maximum, The lower 

limit is given as the reorder level plus a multiple 

of the total demands during a production order cycle 

plus production leadtime period. 

Add the requisition objective of central warehouse 

to total for branches to give product requisition ob-= 

jective. 

Compute the ratio of central warehouse optimal item 

requisition objective to total direct demands of the 

item. This expresses the requisition objective in 

weeks of demands rather than in units, In the main 

program, the requisition objective is updated by mule 

tiplying this factor by the current demand rate of an 

item to give the current level in units again, 

Round off the product requisition objective to the 

next higher unit and compute the product reorder level 

as difference between the product requisition objec-— 

tive and the economic order quantity for the product. 

Sum the product totals to give system totals for: 

1, SD = System Direct Demands Rate: Ail products, 
all items, all locations. 

2.  SDB=- System Branch Demands Rate: Ali preducis, 
all items, all branches except the central 
warehouse, 


on SGP= System Goal for Branches: The System 
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SG PP= 


SDLB- 


SDL- 


SRP= 


SORD= 


SQ= 


SQP=} 


SR= 


Branch demand rate times the branch goal 
period or the number or product units re- 
quired in a specified period (goal). 

System Goel, Plant (Central Warehouse), 

The system direct demand rate times the 
plant goal period or the number of product 
units required in a specified period (goal). 
System Demands in a production Leadtine 

for Branches only. 

System Demands in a production Jeadtime 

all branches and warehouse, 

System Requisition Objective for Branches 
only. 

System Order Rate 

System Economic Order Quantity 

System Economic Order Quantity for Branches 
only. 


System Requisition Objective 


RETURN card = returns the program to the main program 


to commence simulation run on the product and related 


items for which the base rules have just been computed. 


Print Summary Section 


Management Data 


OG = 
V-2 
V=-3 


V=4 


System Goals sum of branch and plant goals, 
Weeks of supply in the branches, 
Weeks of supply in system, 


Difference between number of weeks of 
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oy V-5 
6. Wl - 
7. W2-= 
8. V8 ra 
Ds V9 id 
10. V10= 


supply in the branches under the SRI rules 
and under the alternate rules. 

Difference between the number of weeks of 
supply in the system under the SRI rules 
and under the alternate rules, 

Overall System Branch Goal in weeks, 
Overall System Goal in weeks, 

System Requisition Objective less System 
Demands in a production leadtime. This is 
the approximate maximum system on-hand 
quantity. 

Approximate average on-hand quantity. 
System Requisition Objective for Branches 
less system Demands in production lead= 
time for Branches, This is the approximate 


maximum on-hand quantity in the branches, 


PRINT statement for Base Rules Summary 


FORMAT for PRINT statement of card 815, 


FORMAT statements for the PRINT statements used for 


checkout purposes. 


END (of SUBROUTINE) card, 


SUBROUTINE PLACE BRANCH ORDER (ICON) 


Definition of Subroutine 


Dimension and place in common those subscripted 


variables which are used in other subrcutines of the 


program, 


Note that all subscripted variables used in 
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this subroutine are in common. 

Test to ensure that the quantity to be ordered is 

positive, If quantity is zero or negative no action 

is taken to store the "order" value and progran 
continues. 

Test the branch order quantity against the warshouss 

on-hand balance to see if order can be filled, If 

not, program returns to the mainline, However, if 
the order can be filled the warehouse on-hand balance 
and stock-position are reduced accordingly and the 
on-order tally and stock=-position of the branch is 
increased accordingly. 

An ordere-count tally is kept when using the SRI rules, 

If it is desired to main this tally for other rules, 

change the parameter ICON and read in or set the 

parameter equal to l. 

The following cost and count tallies are updatsd for 

each branch order, 

a. The cost of freight. WJ(J) is the weight per 
unit of the J-th item (lds. per urit); FRK(K) is 
the freight cost per hundred weight for the K-th 
Branch and, of course, the ,O1 is a conversicn 
factor from lbs, to cwts. 

b, The cost of processing (receiving). 

ec, Broken pallet count. This count is expressed in 
the number of units by which the order quantity 


exceeds the full pallet load. UPCJ(J) is the 
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860 


861 


862 
863-881 


number units per case and CPPJ(J) is the number 
cases per pallets, Their product is the number 
of units per pallets, 

d. Broken case count, This count is also expressed 
in units; it is the difference between the order 
quantity and the nearest full case quantity. 

e. Cost of picking and loading. This is the cost 
incident to issuing in non-pallet load or the 
cost of breaking into a full pallet load and 
issuing broken lots, CPFP is the unit cost of 
breaking the full pallet. 

After updating the necessary cost and count tallies 

we are ready to place the branch order and its de- 

livery date into the Branch Order Queue, Here we are 
simulating the delay between the actual placement of 
an order and the availability of the items at the 

Branch, The subroutine parameter, N, is set to 1 to 

route us to the appropriate section of the queuing 

routine. 


End of Subroutine card, 


SUBROUTINE PREDICT 


Definition of Subroutine 

Dimension and place in common those subscripted 
variables which are used in other subroutines of 
the program, No local subscripted variables are 


used in this subroutine. 
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882-383 


884-887 


888 


895-897 


898-899 


The predicted weekly demands are basically a nine 
week moving average, However, since the derand input 
data is monthly demand, a complicated indexing is 
required to select the appropriate month in which the 
most recent nine weeks have occurred and the weighting 
factor given to each. Each monthly demand rate 
selected is scaled down to a weekly rate which is 
then weighted by the number of times the monthiy rate 
is included in the past nine week period, Jl and J2 
are variable indexing limits depending on time, The 
constant values 44 and 52 are chosen so that the 
first weeks prediction will be based on the 11th and 
12th month's demand, (See card 137) 

Initialize total product and total item demand rate 
and set up do-loops on the items and locations, 

Test to see if the branch carries the item. If the 
branch does not carry the item no prediction of 
demands is necessary. The computer will not go 
through the prediction routine for this branch thereby 
saving on computer time, 

If the value of L as computer on card 891 is less 
thai or equal to 12 we test to see if it is zero, if 
it is we add 123; if the value lies between 1 and 12 
inclusively we use that month's demand rate scaled 
down to a weekly value. Nine such weekly demand 
rates are accumulated in the variable X. 


The item=location demand rate is averaged over th= 
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900 


901 


902 


903 


previous nine weeks and the item demand rate is the 
sum of these over all branches, 

Continue Statement for do-loop on branch, If the 
branch does not carry the item, the program proceeds 
to the next branch without waste of time going through 
the prediction routine for that branch, 

Sum the demands for each item to give product demand 
rate, Item demand rates must be converted from units 
to product-equivalent units. 


End card for subroutine, 


SUBROUTINE RDTP (Read Tape) 


Definition of Subroutine, All the input data for 

the product, items, and branches are read into the 

program by this subroutine, The control parameters 

are read by a separate read statement at the begin- 

ning of the program, The following data is read from 

the cards: 

Product Card 

1. The product (Base) identification number 

2, onelf life’ in months 

3, Production leadtime in days 

4. Batch size increment in product units 

5. Batch size mean in product units 

6. Branch Goal (Requisition Level) in days for the 
Alternate rules. 


7, Warehouse Goal (Requisition Level) in days for 
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904~906 


907 


the Alternate rules. 
Item (Size) Card 
l. Size code or item identifasation number 
<. ‘tem conversion factor, i.s., number of product= 
equivalents per unit of the item (¢.2., gellons 
per quart). 
3. Units of the item per case 
4. Cases per pallet 
5 Unit production cost in cents 
6, Unit selling price in cents 
7, Weight per unit (peurds ver unit) times tes 
Branch Card 
1. Branch Code Number 
2. On-hand balance for éach item 
3. Inventory positicn for each item 
4. Twelve monthly demand rates 
Dimension and place in common storage all subscripted 
variables which are used in other subroutines of the 
program, 
READ statement for all thrss types of input cards. 
The input card all have 17 fields and use a comnon 
format statement, The Gata are first placed into a 
temporary storage locatiun (AO(K)); after program 
has returned to the main propram, these data are 
transferred to their permansrt location. This de- 
vice was utilized to standardize the reading of the 


input data. Originally all the product, item, and 
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908 


909 


branch data were pre-recorded on magnetic tapes 

Since the volume of data was too large to read into 
the program efficiently from cards. This procedure 
could be used again by altering this READ statement 
to a READ INFUT TAPE statement, 

Test to see which type of card has just been read. 
AQ(1) is 1 for a product card, 2 for an item card, 
and 3 for a branch card, If a product card has just 
been read the AO(1) flag indicates that it is now 
time to simulate for the previous product, 

If an item card has just been read, a positive value 
in AQ(1) will send the computer back to the main 
program where the subsequent branch cards pertaining 
to the item will be read, 

If a branch card has just been read, a positive value 
in AO(1) will again send the computer back to the 
main program where the program will check for addi-= 
tional items under the same product and for branch 
data under the item, 

After reading the last branch card of the last item 

f the last product, the program again calls the RDRT 
subroutine which attempts to read the next card, Since 
the last input card is a blank; a zero will be read 
into the AO(1) location, The test made on this card 
will cause the value to be reset to a -1. The com-= 
puter then returns to the main program where the 1EOF 


is set to 2 signifying that the last product is ready 
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for simulation and to stop the program after thé 


Simulation is completed, 


910 Return Card 
911 Format for READ statement 
912 End of subroutine card 


SUBROUTINE GENDEM (Generate Demands) 


913 Definition of Subroutine 


‘ 


a 


914-932 Dimension and place in common storagé al. subscr pied 


variables which are used in cther subroutines of the 


1633) 


progran, 
933 An index number is computed +o determine which monthly 


demand rate will be used to compute the actual weexly 


h 


demand rate, The scheduls works out as follows. 


T L Frequency of 
(Time Frame) (Index for Month) Selection 
2. 35450 as : 
6,7;859 2 a 
LOM 12 13 3 PA 
BAP Ake llom ste. Sas: 4 E 
19,20,21,22 5 dy 
PS PIRES LS Plo 6 4, 
Pal Ee oe PAS ho 6 eee a 7 5 
Bic. etc, etc, 
934-935 Set up do=loops on items and lovatiors, 
936 Test to see if branch carries the iter, 
oat Compute the actual weekly demand rats Sy scaiing the 
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938-939 


940 
941-945 


946-963 


964 


monthly demand rate down to a weeKly rate. 
Continue card for do=loops and tnd of subroutine 


card Ry 


SUBROUTINE QUEVE B (N) (Queve Branch Orders) 


Definition of subroutine 

Dimension Local Subscripted Variables, These local 
variables are temporary storage locations for the 
order quantities, delivery dates, item and branch 
incic:s, Once computed they must be stored until 

the proper time to simniate delivery at the branch, 
One hundred and fifty storage locations have been 
reserved for these Branch interim orders, 

Dimension and place in common storage all subscripted 
variables which are used in other subroutines of the 
program, 

The subroutine is divided into three sections; the 
first places the branch order in the queue, the 
second removes it from the queve, and the third sece= 
tion initializes the queve. This computed GO TO 
statement directs the program to cne of three section 
of the subroutine, If the parameter value is 1, the 
subroutine places the order quantity, delivery date, 
and indices in the Queue; if the parameter is 2, the 
subroutine takes the orders out of Queue and simulates 
receipt of order at the branch; if the parameter value 


is 3, the subroutine is initialized in preparation for 
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965 


966 


967-970 


ey 
972-973 


an 7 19 


980-982 


983-986 


Boy 


the next product, 

The I index counts the number of orders in the Queue 
and also identifies the order, 

Test to see if the queue is filled, If it is, this 
information is printed out -;o advise analyst but 
program continues wnich in effect zeros out the order, 
The order, delivery, date, item, and location induces 
are transferred to Queue. 

RETURN card 

FRINT statement for Queue overflow condition and 
RETURN card, 

Test to see if there are any orders in the Queue, 

If there are no orders (I=-0) then the order variable 
(ORD) is zeroed and program returns to the main pro- 
gram, If there are orders in the queve, a do-Loop 
is set up to screen through all the delivery dates 
to see if it is time to remove any orders from the 
queue and to simulate receipt at the branch, 
Transfer the order quantity, item, and location 
indices to main program locations, 

Transfer the last order to quexze to the position 
just vacated by the order whose delivery date was 
due, This is done in order to be able to place the 
next branch order in the last queue position and to 
conserve running time and to minimize the number of 
storage locations for the queue. 


Since the order data of the last queue position has 
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988 


989 


990-991 


oS 


993 


994 


Do e2299 
1000-1017 


1018 


1019-1023 


now been transferred to the vacated slot, the 

order delivery date (QODD) of the last position is 
erased, 

The total number of orders in the queue is reduced 
to account for the order "delivered", 

RETURN card 

This is the initialization section of the queue, 
Format statement for Frint statement used to advise 
of queue overflows. 


End of Subroutine Card 


SUBROUTINE QUEUE P (M) (Queue Production Orders) 


Definition of Subroutine 

Dimension local subscripted variables, 

Dimension and place in common storage those sub-= 
scripted variables which are used in other subroutines 
of the progran, 

The subroutine is divided into three sections; the 
first places the production order data in the queue, 
the second removes it from the queue, and the third 
section initializes the queue, This computed GO TO 
statement directs the program to the appropriate 
section as indicated by the parameter value, 

Use a do=loop to screen the delivery date locations 
for an empty slot to store next production order, If 
all locations are filled, the computer will print 


"Queue overflow" statement and return to the main 
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1024~1029 


1030-1034 


1035-1038 


1039-1040 


1041-1043 


1044-1045 


1046-1048 


program; the order is, in effect, cancelled, 
Transfer the production order data to the queue, 
Allocations to sizes and the allocations to the ware 
house for each item are stored. The machine then 
returns to the main program, 

This is the second section cf the subroutine which 
Simulates the delivery of the production orders, 
Here a do-loop is used to screen the delivery dates 
in the queue to determine if any are due, If there 
are no orders in the gqusue the order is zeroed and 
the computer returns to the main progran, 

The production order data is transferred from the 
queue locations to main program locations, 

The delivery date location of the production crder 
removed from the queue is zeroed, Since there are 
only delivery dates there is no need to maintain the 
orders in consecutive locations as in the Queus B 
routine. 

Initialize the production order delivery dates 
locations. 

Format for Print Statement which advised of queue 


overflows end End-of-Subroutine card, 


MISCELLANEOUS FUNCTIONS 


Definition of RNDNF Function. This function reunds 
the input quantity to nearest lot size, ¢.g., round= 


ing an order of so many units to the nearest case lot 
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1049-1051 


1052-1054 


1055 


288-290 
291-292 


290-200 


quantity, 

Definition of RNDCF Function, This functicn rounds 
the input quantity to the next higher lot size, 
Definition of RNDLF Function. This function rounds 
the input quantity to the next lower lot size. 


End of program card 


ALTERNATE LECISION RULES 


Branch Ordering Procedure heading 

Determine whether it is time to conduct a review of 
the stock records, If it is time, X will be zero, 

If it is not review time then the program jumps to 

the next section, 

Set up a do=loop on items to screen the reorder 

leveis against the inventory positions for each item 
in the warehouse. If the reorder level is less than 
or equal to the inventory position it is assumed that 
the warehouse's stock will be sufficient to meet all 
branch needs and we then test the branch reorder 

level against their inventory position for each item, 
If the branch inventory position is below or equal to 
the reorder level] we compute the order quantity as 

the difference between the branch requisition and 
branch inventory position. If the warehouse inventory 
position had been less than its reorder level indicat-= 
ing that it was low on stock of the item then the 


branch order is reduced by one half, 
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302 


303 


304 
ber er 
308 
poy 312 


313-318 


The delivery date is computed by determining the 

last scheduled delivery date then adding the time 
peried until the next schedule shipment and the ine 
transit time, 

The value of the inventory intransit corresponding 

to this order is added to the previouz total, 

Call the P’ace Branch Order subroutine to simulate 

the filling of the branch order at the central waréq= 
house and updating of the stock records, The PBO sub- 
routine calls the Queve B routine to place the Branch 
order in the queue until the desired delivery date 
arrives, 

CONTINUE card 

Production Ordering Frocedures heading 

Initialize the Production Order Batch Size variable, 
Compute the production quantity as the sum of the 
difference between the requisition level and invene= 
tory position at the central warehouse, Production 
orders are based solely on the item needs at the ware-= 
house under these rules, 

Test to see if production order quantity is positive, 
If it is the value is transferred to the "ord" 

storage location and then the production order quantity 
is rounded to the next higher batch size, To ensure 
that a minimum economic order quantity is produced, 
the crder quantity is taken as the maximum of the 
rounded order quantity and a predetermined mininum 
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A. 


B, 


CHAPTER IV 
OUTPUT REPORTS 


1, Control and branch card input data 


Selectors 


l. 


3. 


4. 


De 


Costs 


1. 


ICON1L. Rule Selector: A value of 1 indicates that the SRI 
decision rules are to be used to compute the requisition 
levels, reorder levels, and economic order quantities; a 
value of 2 indicates that the alternate set of decision rules 
are to be used, 

ICON2. Not used in the present form of the simulator, 

ICON3. Stockout Policy Selector: A value of 1 indicates 
that demands are partially filled, i.e., issues are made 
until a zero balance is reached; a value of 2 indicates that 
demands which can not be fully satisfied are cancelled. 
ICON4,Z. Round to Cartons Selector: A value of 1 indicates 
that Branch interim orders are not rounded to case lots; a 
value of 2 indicates that the order quantity is rounded to 
the next higher full case lot. 

ICON5. Round to Pallet Selector: A value of 1 indicates 
that Branch interim orders are not rounded to pallet lots; 

a value of 2 indicates that the order quantity will be rounded 
to the next higher full pallet lot if it is economical to do 
so, To make this decision the holding cost of the additional 
units is compared to the additional costs incurred due to 


breaking and handling less than full pallet loads, 


Opportunity cost (AZR): The cost of inventory investment in 
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36 


De 


cents per dollar per year. 

Physical Storage costs (APB and APP): The cost attribu- 

table to the physical handling and sheltering of the inventory 
is given in cents per product-unit per year, A different cost 
rate may be assigned to the central warehouse than that which 
is assigned to the branches, 

Order Preparing cost (CPPO): The administrative cost of 
preparing and processing a production order in cents per 
order. 

Picking and Loading (CPL): The additional cost of the physical 
handling of broken cases and pallets is given in cents per 
item-unit handled in a broken case or broken pallet lot, 

Note that there is no additional cost for handling full pallet 
lots since this is taken as the standard handling procedure, 
Cost of Receiving (CR): The costs incurred in the actual 
receiving, handling-in, and storing of items in cents per 


hundred weight. 


C. Multipliers 


1. 


36 


Demand Forecast: This multiplier may be used to increase or 
decrease the predicted demands (see card 898) and thereby 
study the effects of changes in future demand patterns, 

Lead Time: These multipliers are used to study the effects 
of increases or decreases on the branch interim reorder lead= 
times and on the production=-run component of the production 
order leadtime, (See card 76) 

Goals: Multipliers used to vary the branch and central 


warehouse requisitioning goals, (See cards 74 and 75) 
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he 


Standard Cost: Multiplier used to vary the standard unit 


production cost, (See card 102) 


D. Leadtime Components 


i. 


De 


Reaction (LR): The time between the reaching of a reorder 
level and the recognition of that fact and the initiation of 
action to place an order, The input is in days, 

Communcation (LC): The time between reaching a decision to 
place a production order and the receipt of that order by the 
production unit. The input is in days. 

In Transit (LS): The shipping time between the central ware-= 
house and a branch. In the simulator it is assumed for sim-= 
plicity that the shipping time is the same to all branches, 
No provision is made for emergency shipments. The input is 


in days. 


Rules Computation Parameters 


AZR, APB, APP, CR: See section B above, 

IR, LC, LS: See section D above. 

Al: The cost of setting up a production run in dollars per 
set up. 

A3: Not used in the present form of the simulator, 

MU: The system inventory constraint factor, This multiplier 
is used to reduce the "optimal" requisition objective of each 
item-location by some multiple of the demands in a branch 
production leadtime. (See cards 707 and 782), It is also 
applied to reduce the economic order quantity (see card 665). 
By varying the value of MU the system average holding cost or 


system average on hand quantity can be set to any pre-determined 


qd 


o~ 
—_ - = 
‘se <4 
- = = > —- 
— ee _—S —- 
1 LL. a — « — - 
- - ; © lap, 
> —- «& — —- =) es © 
== © 
——_ | —_ 
a - — * @t o' 
> “as 
-— aaa _ »_ = —_- «= aap @& 
—_ 
—_ . = , Ss © 
— ——- > = ——> ss 
e — 
- —=—=—=- 
> =a ——— eae —_— = - 





valve. Setting MU equal to zero is equivalent to placing no 
constraints on the system, 

K6: Multiplier used to increase or decrease the central 
warehouse item reorder level in the RXCOMP subroutine, (See 
card 755). 

K5: Multiplier used to increase or decrease the system 
item-reorder level in the RXCOMP subroutine, (See card 751), 
K4: Multiplier used to increase or decrease the item-location 
minimum requisition objective, The minimum level is taken as 
the reorder level plus the demands during the production cycle 
plus production leadtime period, The multiplier increases 
only the later component of the minimum requisition cbhjective, 
Hence, doubling the multiplier will not necessarily double the 
minimum requisition level, This multiplier applies to the 
central warehouse's minimum level as well as to the other 
branches, (See cards 720 and 784 in RXCOMP subroutine), It 
should be noted that this multiplier is used as a dummy vari- 
able within the program and is not an input quantity. Frior 
to computing the "optimal" requisition objective level for the 
central warehouse, the value of the multiplier FK3 is put into 
FK4 (card 768), and prior to computing the level for the other 
branches the value of FK1 is transferred to FK4 (card 680). 

It is not known why the variables FKL and FK3 were not used 
directly in computing the minimum requisition objectives, 

K3: This multiplier is used in the main program (card 145) 

to update the "optimal" central warehouse requisition level. 


It is used to increase or decrease the lower limit on the 
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ll. 


ue 


13. 


14. 


requisition level, See K4 for further information, 

K2; Multipiier used to increase or decrease the item-location 
reorder level for all locations except the central warehouse 
(see card 640, 719, and 151). By varying this multiplier the 
lower limit on the requisition objective is also changed since 
the minimum requisition level depends on the reorder level, 
Kl; This multiplier is used in the main program (card 153) 

to update the "optimal" item-location (except for central 
warehouse) requisition objective, It is used to increase or 
decrease the lower limit on the requisition level, See K4 

for further information, 

CM2: Cost multiplier used to increase or decrease the pro- 
duction setup cost (card 665). 

VLPl:; The rel-variance (ratio of variance to the square of 
the mean of the prodiction=run component of leadtime), It is 

a dimensionless quantity. 


VLP3: Not used in the present form of the simulator, 


F. Branch Data 


1. 


36 


Ae 


Codez: The activity accounting numbers assigned to the 
branches, Code 1 is assigned to the central warehouse. 

Cost of Shipment: The freight rate in cents per hundred 
weight, The freight rate at the central warehouse is given 
the negative value of the receiving cost in order that the 
unit cost at the warehouse will be the unit production cost, 
Scheduled Shipping Frequency: The period between shipments 


to each branch. 


Branch Inventory Constraints Multipliers: If it is desired 
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to limit the individual branch total average inventory tc 
some specified dollar value or on-hand quantity, vary this 
parameter on a trial and error basis until the proper setting 
is obtained. Each branch can be given a different constraint 


factor if it is so desired, 


2. The item, product, and system summary report 

The average inventory values, sales, costs, and stock effectiveness 
data are summarized here by item, product, and system, Complete branch 
data is listed for each item. The "Holding Cost" given in the report 
includes the opportunity cost factor as well as the physical holding cost, 
The "Order Count" is the number of interim branch orders for the item. 
The "Estd Cost" of stockouts is the profit loss due to lost sales, The 
number of broken cartons and pallets is given in item units rather than 
in the number of cartons or psllets actually opened to fill an order, 


The other headings on the report are self-explanatory. 


3, The branch summary report 

The inventory values, cost data, and stock effectiveness measures 
are summarized by branch for all items and products, 

All the data given in this report under the "Inventory Values" 
heading including the stock turn ratio are expressed in terms of dollars. 
The "Order Count" is the number of scheduled deliveries to the branches, 
The "Preparation and Processing" cost is the cost of making the scheduled 
deliveries to the branches, It is the number of scheduled deliveries per 
year times the cost of preparing and processing an order (CPPO), The 


"Picking and Loading" cost is the sum of the costs incident to handling 





of broken cartons and broken pallet lots and the cost of handling the 


full pallet lot portion of the orders, The "Unloading Processing" cost 


is the cost of receiving and storing materials at the branch, The 


other headings on the report are self-explanatory, 


3. 


he 


De 


10, 


ll. 


4. The base rules summary report 
D: The expected demands per week for the system, 
RPR/D: The requisition objective in weeks for all branches except 
the central warehouse under the SRI rules, 
R/D: The requisition objective in weeks for all branches including 
the central warehouse under the SRI rules. 
APR/D and A/D: Same as 2 and 3 above except these are for the 
alternate policy rules, 
(RPR-APR)/D and (R-A)/D: The difference between the requisition 
objective requirements under the SRI and the alternate policy rulés. 
DL( PR): The "pipeline" or leadtime quantity in product units for 
all branches except the central warehouse, 
(DL): The "pipeline" or leadtime quantity in product units including 
the central warehouse. 
(ON-HND=MX): The average maximum on-hand level for all branches 
including the central warehouse. 
(ON-HND=MX)PR:e Same as 8 above but does not include the central 
warehouse. 
(ON-HND=AVE): The expected on-hand inventory level for all branches 
including the central warehouse, 


(ON-HND-AVE) PE; Same as 9 except the central warehouse is not 


included. 
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ORD: The average number of production orders per week. 

Branch Picking Cost: The system-wide cost of picking and loading 
for both branch interim orders and production allocationg, 
Production Cost: The total set-up costs. 

Shortage Cost: The loss of profit due to nonavailability of items 
to meet demands as they occurred, 

Holding Cost: The physical holding costs incurred at all branches 
including the central warehouse, 

Average Value in Transit: The dollar valve of the average number 


of units in transit between the central warehouse and the branches, 
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CONTROL CARD INPUTS -- RUN NUMBER 5 


SELECTORS -- RULE L 
STOCKOUT 1 

ROUND TO CASES (ad 

ROUND TO PALLETS 2 

COSTS -- OPPORTUNITY «2065 
; , PHYSICAL STORAGE - BRANCHES 024 

- WAREHOUSE 043 

ORDER PREPARING — 075. 


PICKING AND LOADING = BROKEN PALLETS ,.07 
- BROKEN CASES 020 
- FULL PALLETS 012 


RECEIVING 15 

MULTIPLIERS -- DEMAND FORECAST 1.00 
LEADTIME - WHS TO BRANCH 1.00 

- PLANT TO WHS 1.00 

GOALS - BRANCH 1,00 

- WAREHOUSE 1.00 

STANDARD CO&T 1.00 

LEADTIME -- REACTION 46 
COMMUNICATION l. 

IN TRANSIT ll. 


PASE RULES COMPUTATION 


AZR 0205 APB ~240 APP 430 CR «150 
LR 4.000 LC 1.000 LS 11.000 Al 400 .000 
AS ~000 MU 000 K6 1.000 K5 1.000 
K1 1.000 K2 1.000 K3 1.000 C2 1.000 


VLP1 ~000 VLPS 1.000 





BRANCH DATA 


CODE 


CS 


MUK 


1. 45. 57. 62. 63. 64. 68. 
~.15 83 042 oD 059 oor or 
0. La is as ls » dee le Le 


«0.0 «0 .0 20 .0 70720 50 


88. 106. 107. 1ll. 122. 127. 136. 138. 150. 175. 
087 .56 34 226 .23 4.32 0435 049 .32 235 147 
Ze bie Ce le Ve Se oe Se Ze ie Ze I i. 1. Aes 


Pas) 20 ALS) sO AS, eC «0 0 A®) 0 <0 AG; P18, 0 ae) 


180. 





Gara l. 


CHAPTER V 


CONTROL CARD FORMATS 


This card contains the run number, the control selectors 


settings and the system cost factors, 


Ge 


Ce 


Be 


Columns 1=5, The letters "IDENT" are placed in these columns 
to identify the run number, If desired these columns and/or 
any columns on subsequent cards used for alphabetical identi- 
fication may be left blank, 

Columns 6-8, These are usually left blank but may be utilized 
for any non-input data as the letters "IDENT" in columns 1-5, 
Columns 9-14. The run number designator is entered here, It 
may be an alphanumeric designator, 

Columns 15-24. The letters "SELECTORS" are placed in the 
columns, right adjusted to leave column 15 blank, 

Columns 25=26, 27-28, 29=30, 31=32, 33=34. The values for 
ICON1, ICON2, ICON3, ICON4, and ICON5 (i.e., the number 1 or 
2) are placed in these columns, right adjusted, 

Columns 35-40. The letters "COSTS" are inserted here, also 
right adjusted. 

Columns 41-44. The opportunity cost (AZR) is inputed in 

mils, right adjusted. A value of 205, for example, indicates 
an opportunity cost of 20.5 cents per dollar of investment per 
year. 

Columns 45-58, The physical holding cost for the branches 
(APB) is inputed in cents, right adjusted, A value of 24, for 


example, indicates a holding cost of 24 cents per procuct unit 
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Card 


Me 


De 


Q6 


per year, 

Columns 49-52, The physical holding cost for the central 
warehouse (APP) is inputed in cents, right adjusted, A valus 
of 43, for example, indicates a holding cost at the warehouse 
of 43 cents per product unit per year, 

Columns 53-57, The cost of preparing and processing an order 
(CPPO) is placed here in cents, right adjusted, The value of 
75, for example, indicates a cost of 75 cents per order, 
Columns 58-62, The cost of handling broken pallets (CBPC) is 
given here in cents, right adjusted, The value of 7, for 
exampie, indicates a cost of 7 cents per item unit handled, 
Columns 63-67, Same as k above except for broken cas#s (CBCC), 
Columns 68-72, The cost of receiving (CR) in cents per 
huncred weight, right adjusted, 

Columns 73-77, The cost of handling (or picking) a full 
pallet lot (CPFP) given in cents per pallet, right adjusted. 


This card contains the management multipliers, the reaction, 


communication, and shipping components of the production leadtime and 


the cost of a preduction set-up, 


Be 


Columns l=-ll. The letters "MULTIPLIERS" is inserted here for 
identification of the card and the subsequent input data, 
Columns 12-15, 16-19, 20223, 24=27, 28-31, 32=35, The values 
for the demand forecast (FMDF), the warehouse to branch (FMLWB), 
and the plant to warehouse leadtime (FMLP) multipliers, ths 
branch and plant goal multipliers (FMPG, FMBG), and the stan- 
dard cost multiplier FMCS are placed in these columns in the 


given order, The multipliers are specified to two decimal 
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Card 


Ee 


3. 


places within the program; hence, a value of 100 indicates 

a multiplier of 1.00. A setting of 1.CO for the multipliers 
is equivalent to not using them, They should never be set at 
zero. 

Columns 36-45, Leave blank 

Columns 46-55, The words "LEAD TIMES" are inserted here, 
Columns 56-58, 59-61. The value of the reaction (FLR) and 
communication (FLC) components of leadtime are given here in 
days, right adjusted, 

Columns 62-64, Leave blank 

Columns 65-67, The value of the shipping or in-transit 
component of leadtime (FLS) is given here in days, right ad- 
justed, 

Columns 68-72, The letters "Al" indicating the set-up cost 
are placed here, right adjusted, 

Columns 73-75. The value of the set-up cost is inputed here 
in dollars per set up, right adjusted. 


This card contains the input data for the system capacity 


constraint (FMU), the reorder level multipliers (FK6, FK5), the 


minimum requisition objective multipliers (FK1, FK2, FK3), which are 


referred to ae "floors" in 1. » the set up cost multiplier (CM2), 


and the rel-variance for the production-run component of leadtime 


(VER). 


Be 


Columns 1=3, The letters "MU", left adjusted, are listed here 


for identification only. 
Columns 4-8 The value of the system capacity constraint (FMU) 


is given here to three decimal places (i.e., a value of 8300 
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Co 


CG. 


Ge 


Card 4. 


(Note: 


is read in as 8,300), 

Columns 9=24, The letters "TRIGGER LEVELS" are entered here. 
right adjusted, for identification only. 

Columns 25=27, 28-30, The values of FK6 and FK5 are entered 
here to one decimal place, (i.e., a value of 14 is read as 
1.4), They should never be set to zero, 

Columns 31-38, The letters "FLOORS" are placed here, right 
adjusted, for identification only. 

Columns 39-41, 42-44, 45-47, The values of FKL, FK2, FK3 

are entered in these fields in their respective order, The 
multipliers are given to one decimal place (i.e., a value of 
O60 is read as 6.0). If these multipliers are set at 1.0 it 
is equivalent to notutilizing them, The multiplier FK2 
should never be set to zero, 

Columns 48-59. The letters "VARIATIONS", right adjusted, are 
placed here for identification only. 

Columns 60-63, The value of the set-up cost multipliers (CM2) 
igs entered here to two decimal places, (i.e., a value of 404 
is read as 4.04). 

Columns 64-71, The value of the rel-variance of the produc= 
tion-run component of the production leadtime is entered here 
to two dscimal places, right adjusted. A value of 150 is 
read as 1.50. 

Branch Data heading card to identify subsequent input cards. 


The READ format statement is designed to skip a card at this 


point and would mis=read all the subsequent cards if this particular 


one is omitted) 
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Card 


Card 


Card 


Card 


Card 


Ds Branch Code Card, Three columns are used for each branch 
code with the code numbers right adjusted, Any set of code numbers 
(each not exceeding three digits) may be used, There is a limita-— 
tion of 25 branch code numbers. 
6. Freight Rate Card, The cost of shipping,FRK(K), in cents 
per hundred weight is entered for each branch in same order as the 
branch codes are entered, Three columns are allowed for each freight 
rate cost. The limitation is also 25 entries, with the numbers rizhi 
adjusted, 
ia Shipping Schedule Card, The scheduled shipping frequency for 
each branch, FSK(K), is listed on this card, (i.e., the number of 
weeks between shipments from the central warehouse to the individual 
branches, Three columns are allowed for each field, The values 
are given as integers and are right adjusted. 
8. The Branch Capacity Constraints, The individual branch con~- 
straint multipliers (FMUK(K)) are entered here to one decimal places. 
Three columns are allowed for each entry and numbers are right 
adjusted. A value setting of zero is equivalent to placing no 
constraints on the branches. 
a Product Data Card 
Be Column 1, The number 1] is always placed here to indicats 

that this is a product data card. 
ler Columns 2-5, The base identification number of the product 

is entered here. 
on Columns 6-20, Leave blank 


d. Columns 21=-25. Enter the shelf-life in months, right adjusted, 


Ce Columns 26-50, Leave blank 


8&9 


Card 


nO’. 


Ge 


Columas 51-55, Enter the length of a production run (FLP) in 
days, right adjusted, 

Columns 56-60, Enter the batch size increment (BSI) in 
product units, right adjusted, 

Columns 61-65. Leave blank, 

Columns 66-70, Enter the minimum economic batch size (BSM) 

in product units, right adjusted, 

Columns 71-75, The Branch Goal (BG) is given in days, right 
adjusted. 

Columns 75-80, The Plant Goal (PG) value is given in days, 
right adjusted. 

Item Data Card 

Column 1, The number 2 is always placed here to indicate tha* 
this is an item data card, 

Columns 6-9, The product base identification number is entered 
here, right adjusted. 

Columns 10-11, The item identification number is entered here, 
right adjusted. 

Columns 15=20, The item=-unit to product-units conversion 
factor (or size factor) is entered here to four decimal places, 
The value of 11057 is read as 1.1057. 

Columns 51-55, The number of units per case (UPCJ(J)) of the 
item is entered here, right adjusted. 

Columns 56-60, The number of cases per pallet for the item 
(CPPJ(J)) is entered here, right adjusted. 

Columns 61-65. The unit production cost (CJ(J)) of the item 


in cents is entered here, right adjusted, 
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Card 


C. 


Columns 66-70, The sales price of the item per unit (PJ\J)) 
is entered here, right adjusted, 

Columns 71-75, The weight per unit (Wi(J)) of the item in 
pounds is entered here, right adjusted. 

Notes: All other columns are left blank 

The Item=-Location Data Card 

Column 1, The number 3 is always entered here to indicate 
that this is an item-location data card, 

Columns 12-14. The branch code (BCK(K)) is entered here, 
right adjusted. 

Columns 15=20, The initial on hand balance in item units 
(OHJK(J,K)) is entered here, right adjusted, 

Columns 21-25, ===, 76-80. These are 12 fields of five 
columns each in which the monthly demands are entered in 
item units, The average long run demands based on these 
numbers will be used to make the initial forecast of demands 
and to compute the initial optimal inventory policies, They 
are also used individually to compute the actual weekly 
demands lodged at the branches, 


Notes All other columns are left blank, 
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CHAPTER VI 


PROBABILISTIC DEMAND PATTERNS 


Four subroutines have been designed to utilize probabilistic demand 
patterns for the SRI simulator. These subroutines will provide 12 random 
numbers from the appropriate distribution desired for each of the monthly 
demands. The procedures for scaling the monthly demands to arrive at the 
predicted and actual weekly demands remain the same. 

Normal, exponential, uniform, and Poisson distributions are approxi- 
mated in these subroutines. The distribution parameters are defined 
within the subroutine and can be changed to any desired value by changing 
the appropriate cards, The distributions are truncated to prevent un= 
realistic values of demand from occurring, In the case of the continuous 
functions, the random numbers are rounded to the nearest integer, Even 
with these modifications a close approximation to the real distributions 
has been attained which is more than accurate for the purposes intended. 
The normal distribution is truncated two standard deviations above and 
below the mean, The exponential distribution is truncated at that value 
above which only a certain desired probability (error term) occurs. The 
uniform distribution, of course, is bounded by its upper and lower end 
points which may be set to any desired values. The Poisson distribution 
is terminated at three times its mean, Table 1 gives the results of 
some test runs on these subroutines, 

In order to use the random demand subroutines the program and input 
cards must be modified as follows: 

iL. Input Cards, On each branch data card leave columns 21 through 


80 blank instead of inserting the monthly demand values, 
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4. 


Subroutine Parameter Cards, Make new cards with the desired 

values of the distribution parameters and in the case of the 

exponential subroutine set the desired value for "ERROR", 

i.,e., the probability area to be truncated (a value of .01 to 

.05 is suggested). 

COMMON Cards, All COMMON cards of the main program must be 

listed in the subroutine and similarly all COMMON cards of 

the subroutine must be listed in the main program, 

Call Cards 

a. Initialize the Random Number Generator, [In order to vary 
the random numbers used in different runs the random de- 
mand subroutine may be called a few times at the beginning 
of each run by the following routine: 


ISPIN = ( ANY DESIRED FIX POINT NUMBER) 65A 


CALL UNIFORM (1) 65B 
DO 99 K = 2, ISPIN 65C 
99 CALL UNIFORM (2) 65D 


b, Call for the Monthly Demands, Insert the distribution 
call card after card number 118; for example: 


CALL NORMAL (2) 118A 


The above procedures will introduce only one random demand pattern 


for all products and items since the distribution parameters are fixed 
for each simulation run, If it is desired to vary the parameters for 
different products, items, or locations, this may be accomplished by 
listing the parameters on either the product, item or branch data card. 
Since each distribution has two parameters columns 26-30 and 31-35 of 


the data cards may be used for this purpose. The RDTP subroutine will 
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place these values in the A0(15) and AO(16) locations, 


Products: To vary the parameter values for each product, use the 


following routine: 


ie 


3. 


Items: 


Insert the value of the two parameters in columns 26-30 and 
31-35 in the product data card, right adjusted, These values 
will be read as integers and must be scaled if decimal places 
are desired, 

Place the names of the parameters in COMMON and remove the 
parameter cards from the subroutine. 

Prepare cards for the main program to transfer the parameter 
values from the AO array to their permanent storage locations. 
For example, if using the SUBROUTINE NORMAL, prepare the 


following cards: 


FMEAN = A0(15) 734 
STDDEV = A0(16) 73B 
STDDEV = STDDEV * .02 73C 


Insert these cards after card 73 in the main program, The 
last card is an example of scaling a value to a desired decimal 
figure and may or may not be required, 


To vary the parameter values for different items under the 


same product, use the following procedure: 


1. 


Same as (1) for products except insert parameters values on 
the item data cards instead of the product data cards. 

Same as (2) for products 

Same as (3) for products except number the transfer cards 


98A, 98B, and 98C and insert these cards after card 98 in 


the main program. 


ye 





Items-locations: To vary the parameter values for each location 
which carries the item, use the following procedure; 

ly Same as (1) for product except insert the parameter values on 
the branch data card instead of the product data card, 

aa Same as (2) for products, 

3. Same as (3) for products except number the transfer cards 
116A, 116B, 116C and insert these cards after card 116 in 
the main progran, 

SUBROUTINE UNIFORM, Uniformally distributed random numbers are 


generated on the interval (A,B) by the formula, 
U(A,B) = A + (Be-A)*#RN 


where RN is a U(0,1) random number produced by SUBROUTINE RAND. The 
above formula was derived by equating the integral of the uniform density 
between the lower end point and a variable upper limit to the U(0,1) 
random number generated by RAND and solving for the variable upper limit. 
The variable "FLWPT" is the lower end point in the floating point mode 
and the variable "UPENDPT" is the upper end point, also in the floating 
point mode, These variables may be set at any desired values. 

SUBROUTINE EXPONTL. Exponentially distributed random numbers are 


computed using the formula, 


bé 


zl Log, ( 1./(1. = RN) ) 


where RN is a U(0,1) random number produced by SUBROUTINE RAND. This 
formula was also derived by the integral method, The variable "THETA" 
is the parameter of the distribution and not the mean which is, of 


course, the reciprocal of the parameter, The variable "ERROR", used 
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elsewhere in the subroutine, is the value of the probability area under 
the density curve truncated in order to rule out unrealistic demands. 
SUBROUTINE NORMAL. The normally distributed numbers are computed 
by the well known method of summing 12 U(0,1) random numbers, subtract- 
ing the number six, and multiplying this quantity by the desired standard 
deviation and then adding the desired mean, The variable "FMEAN" and 
"STDDEV" are the mean and standard deviation desired for the distribution. 
Since it is possible to generate negative demands with this distribution 
the standard deviation and mean chosen should be such that the mean less 
twice the standard deviation is greater than zero, On the other hand, 
negative demands could be allowed to simulate materials returned to 
store or items to be exchanged. 
SUBROUTINE POISSON, The discrete analog of the integral method is 
used for this distribution. The probability mass function is summed 
term by term and compared to the U(0,1) random number obtained from the 
RAND subroutine. The Poisson number is taken to be the integer correspond= 
ing to that term which when added to the sum of the previous terms ex= 


ceeds the U(0O,1) random number, The variable "FLAMBDA" is the mean of 


the distribution. 
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30 


35 


SUBROUTINE NORMAL (ICON6) 


COMMON INDX, NX, NRNC, NXIN, NRC, NSTRN 
INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 


GO TO (1,2), ICON6 


ITRAN = 1 
NRNC = 10000 
NRC = O 
NSTRN = O 
INDX = 6 
FMEAN = 50.0 


STDDEV = 25,0 
UPLIMIT = FMEAN + 2. *STDDEV 
FLWPT = FMEAN <2, *STDDEV 


DO 38 J =1,12 


M = 13 +J 
Y = 0,0 
DO 5 1=1,12 


CALL RAND (ITRAN) 

RN = FLOATF(NXIN)/10000. 

Y = Y + RN 

AO(M) = FMEAN + STDDEV*(Y-6. ) 
IF(AO(M) =-FLWPT) 30, 30,35 
AO(M) = FLWPT 

GO TO 38 


IF(AO(M) = UPLIMIT) 38, 36, 36 
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36 AO(M) = UPLIMIT 


38 AO(M) = INTF( AO(M) +.5) 
RETURN 


END 
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SUBROUTINE EXPONTL (ICONG) 


INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 
COMMON INDX, NX, NRNC, NXIN, NRC, NSTRN 


GO TO (1,2), ICON6 


ITRAN = 1 
NRC = 0 
NSTRN = O 
INDX = 6 
NRNC = 10000 
THETA = ,Ol 
ERROR = 0,01 


UPLIMIT = LOGF ( 1./ ERROR ) / THETA 
UPLIMIT = INTF( UPLIMIT ) 
Deroad=leae 


CALL RAND (ITRAN) 


RN = FLOATF (NXIN) / 10000, 

M=13+d 

AO(M) = INTF( LOGF(1./(1. - RN))/THETA +0.5) 
AO(M) = MIN1F ( AO(M), UPLIMIT ) 

RETURN 

END 
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SUBROUTINE POISSON (ICON 6) 


INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 
COMMON INDX,NX,NRNC ,NXIN,NRC ,NSTRN 
GO TO (1,2), ICON6 

ITRAN = 1 

mpX = 6 

NSTRN = 0 

NRC = 0 

NRNC = 10000 

FLAMBDA = 4. 

PLE = 12 

Waetl J = 1,12 

CALL RAND (ITRAN) 

RN = FLOATF (NXIN) / 10000. 
TESTNR =RN ¥*EXPF(FLAMBDA) 
M=13+J 

IF(TESTNR=1.)5,5,8 

AO(M)=0,0 

GO TO 11 

PRODUCT=1.0 

1.0 

DO 9 K=1, IUPLIM 

XNR1 = FLOATF (K) 


PRODUCT = PRODUCT *XNR1 
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Y = FLAMBDA**K/PRODUCT + Y 
IF(TESTNR-Y)10,10,9 

9 CONTINUE 

10 AO(M) =FLOATF( K ) 

11 CONTINUE 
RETURN 


END 
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SUBROUTINE UNIFORM (ICON6) 


COMMON STORAGE 
INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 
COMMON INDX, NX, NRNC, NXIN, NRC, NSTRN 


GO TO (1,2), ICON6 


ITRAN = 1 
NRC = 0 
NSTRN = O 
INDX = 6 
NRNC = 10000 
FLWPT = 50.0 


UPENDPT = 80.0 

DOm5e0 = lye 

CALL RAND (ITRAN) 

RN = FLOATF (NXIN) / 10000, 

M=13+d 

AO(M) = INTF( FLWPT +(UPENDPT - FLWPT)*RN + 0.5) 
RETURN 


END 


103 





10 


20 


SUBROUTINE RAND (ITRAN) 


RANDOM NUMBER GENERATOR 


LOCAL STORAGE 

DIMENSION XI(8) ,NXI(8) 

INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 
COMMON INDX,NX,NRNC ,NXIN,NRC ,NSTRN 


EQUIVALENCE (XI,NXI) ,(11,NT1),(INDX,DX) ,(D1,ND1) 


GO TO (10,20), ITRAN 
ITRAN=2 

¥I(1)=3514 6016 2524 6131 
X¥1I(2)=0337 1363 2712 7740 
XI(3)=1760 3011 0710 3016 
X1(4)=0670 1154 5656 1316 
X1(5)=1506 7663 7414 1566 
X1(6)=2100 7160 3140 2631 
X1(7)=1037 4327 7340 4007 


X1(8)=0 


NRC=NRC+1 
L=INDX+1 
ND1=L 


DX=D1 *7 
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30 


M=IND¥+4+1 

ND1=M 

Ti=Dee7 

N=NT141 

NXI(M)=NXI(L)4+NXI(N) 
X¥I(M)=X1(M) *3777 7777 7777 7777 
NX¥=NXI(M) 


LDA(NX) ,MUF(NRNC) ,STA(NXIN). 


RETURN 


END 
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TAPLE 1 


RESULTS FROM TEST RUNS ON THE RANDOM DEMAND SUBROUTINES 


NORMAL DISTRIBUTION SAMPLE SIZE 500 


MEAN XBAR STDDEV SAMPLE VAR Vs. ve 
25 26,54 25 550.1 ood, 
50 49 Lh 25 615.6 24.8 
75 74.79 25 561.3 2nd 

100 101,71 25 547.7 220% 
125 124534 25 506.99 22.5 
150 149.50 25 511.69 22.6 
EXPONENTIAL DISTRIBUTION SAMPLE SIZE 500 
MEAN XBAR STDDEV SAMPLE VAR Sa 
100 102.6 100 9,657.9 98.30 
20 19.78 20 357 82 18.99 
10 9.89 10 88, 86 9.42 
i; Oy 4 15.16 3.90 
2 1.88 2 3.46 ie 3 
POISSON DISTRIBUTION SAMPLE SIZE 500 
MEAN XBAR STDDEV SAMPLE VAR S. Vv. 
4 4.10 2.00 3.99 2.00 
5 5.06 2e27, se 2.29 
6 5.92 2065 5.10 2.26 
7 7.10 2.64 ees 2.67 
g 7.83 2583 Tool Quad, 
9 9.06 3,00 05h 3.09 
10 9.98 3816 9.74 a .12 
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UNIFORM DISTRIBUTION SAMPLE SIZE 500 


INTERVAL EAN XBAR STDDEV SAMPLE VAR Vs. Ve. 
(0,100) 50 50.91 28.86 823.59 28.69 
(10,90) 50 49.81 23.08 526.92 22.95 
(20,80) 50 49.07 aoe 273.34 16.53 
(30,70) 50 Lae, 11.53 Ile a eze: tee? 
(40,60) 50 49.13 5.77 32.13 5.66 
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Variable 
Name 


A 
AQ 
APB 


Ber 


AZR 

BASE 

BCCJK J,K 
BCK K 
BCP 

BCFCJ J 
BCPCJK J,K 
BCPI 

BPCJK J,XK 
BPIJK 

BSI 

BSM 

BF J,k 
BG 

BPCI 

CBC 

CBCC 

CBCI 

CBCJ 

CBCJK 

CBP 

CBFC 


APPENDIX II 


DICTIONARY OF VARIABLES 


Definition 


Production setup cost 
Input data array 
Branch physical storage cost 


Whs physical storage cost 


Opportunity cost 

Product code number 

Broken case count, item=-location 

Branch code 

Broken ctn and pallets, system 

Broken ctn and pallet count, item 

Broken ctn and pallet count 

Broken ctn and pallets, product 

Broken pallet count, item=-location 

Computational factor 

Batch size increment 

Batch size, minimum 

Branch usage flag 

Branch requisitioning goal 

Broken pallet count, product 

Cost of broken ctns, system 

Cost of picking broken ctns 

Cost of broken ctns, product 

Cost of broken ctns, item 

Cost of broken ctns, item-location 

Cost of broken pallets, system 

Cost of picking broken pallets 

Cost of broken pallets, product 

Cost of broken pallets, item 

Cost of broken pallets, item-loca= 
tion 

Cost of branch orders, system 

Cost of branch orders, branch 

Freight cost, system 

Freight cost, branch 

Holding costs, system 

Holding costs, product 

Holding costs, item 

Opportunity plus physical holding 

Holding costs, item=-location 

Holding cost rate, branch 

Whse unit cost 

Branch unit cost 

Production unit cost 

Cost of losses, system 

Cost of losses, product 


Ie 


Unit 


Dollars/Setup 

None 
Cents/Product unit/ 
Year 

Cents/Product unit/ 
Year 
Cents/Dollar/Year 
Noneé 

Unit 

None 

Unit 

Unit 

Unit 

Product Unit 

Unit 

None 

Froduct.Unit 
Product Unit 

None 

Days (of Supply) 
Unit 

Dollars 
Dollars/Unit 

Dollars 

Doilars 

Dollars 

Dollars 
Dollars/Unit 
Dollars 

Dollars 


Doliars 

Dollars 

Dollars 

Dollars 
Dollars/100 Pounds 
Dollars/Year 
Dollars/Year 
Doliars/Year 
Dollars/Unit/week 
Dollars/Unit/Week 
Doliars/Unit/Week 
Dollars/Unit 
Dollars/Unit 
Dollars/Unit 
Dollars 

Dollars 


a 





Cid 
CLJK 
CLK 
CKZJK 


CP 
CPD 
CPDK 
Cae 


CPPK 


CPK 

CPFP 
CPLK 
CPPJ 
CPFO 


CPIJK 
CR 


CSIJK 


CTOT 
CSZ 
CFL 
DBI 
NBIJJ 


DIJ 
DIJK 
DICHI 


DIFR 
DIPRFR 
DLBI 


DLI 
DLIJ 
DM 


DSCAJ 
DSIJ 

DSIJK 
DWIJK 


DZIJK 


J, 


K 


Cost of losses, item 

Cost of losses, item=-location 

Cost of losses, branch total 

Variance factor 

Order rate multiplier 

Cost of setups, system 

Cost of loss per demand, system 

Cost of loss per demand, branch 

Cost of order preparation and 
process, system 

Cost of order preparation and 
process, branch 

Processing cost, branch 

Cost of picking full pallet, rate 

Cost of picking and loading, branch 

Cartons per pallet 

Order preparation and processing 
cost, rate 

Processing cost, item-location 

Cost of receiving (processing) per 
unit 

Profit (selling=branch cost) also 
cost of shortages for both direct 
shortages at branches and short= 
ages resulting from stockout at 
the central warehouse 

Sum of CPLeeCP, ChyesH, SYS=-TOT 

Cost of losses, item 

Cost of picking, loading, system 

Demand rate, branches less warehouse 

Demand rate, branches less warehouse 


Predicted demana, product 
Predicted demand, item 
Predicted demand, item-location 


Mean holding cost/unit of product 
times system demands(DI) 

Mean demand between orders, product 

Mean demand between orders, item 

Leadtime demands, branches less 
whse, product 

Leadtime demands, branches less 
whse, item 

Leadtime demands, whse, product 

Leadtime demands, whse, item 

Monthly demands 

Probable demand in half of leadtime 

Discrete adjustment interval 

Satisfied demand, item 

Satisfied demand, item-location 

Weekly demands, item-location 

Sum of DZIJK over all products 

Mean of distribution for number of 
demands in prodcn cycle plus 
leadtime 
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Dollars 
Dollars 
Dollars 

None 

None 

Dollars 
Dollars /Unit 
Dollars/Unit 


Dollars 


Dollars 
Dollars 
Dollars/Pallet 
Dollars 
Ctns/Pallet 


Dollars/Order 
Dollars 


Dollars/100 Pounds 
Dollars/Unit 


Dollars 

Dollars 

Dollars 

Product unit/Week 
Unit/week 

Product unit/Week 
Unit/week 

Unit/Week 
Dollars/Product unit 


Unit 
Unit 
Product unit 


Unit 


Product unit 
Unit 
Unit 
Unit 
Week 
Unit 
Unit 
Unit 
Unit 
None 
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F Correction factor to Wilson's EQQ; Product unit 
Sum of all FIJK 

F2 Square of F Product unit squared 

FITEM J Item number, floating point mode None 

FIJK The demands in a production iead= Unit 
time adjusted for the variances 
of the demand and leadtime dis-= 
tributions 

FKL Multiplier, branch requisition None 
objective floor 

FK2 Multiplier, branch reorder level None 
{oor 

FK3 Multiplier, warehouse requisition None 
objective floor 

FKL Multiplier, warehouse reorder level None 
Toor 

FK5 Multiplier, item production trigger None 
level 

FK6 Multiplier, production trigger None 
level, warehouse 

FKIJK Variance factor for demands and None 
production reorder leadtime; 
used to compute RIJK(J,K) for 
branches other than warehouse 

FKZJK Same as FKIJK except applies to None 
branch interim reorder leadtime; 
used to compute RIJK(J,1) 

FL Warehouse production leadtimes; for Week 
central whse LC + LP + LRI/2 

FLC Communication leadtime Week 

FLIK Branch production leadtime; Sum of Week 
IRI/2, 10 , LE, ISK, bWy 2 

FLDI Maximum BOQ due to shelf life Product unit 
limitations 

FLP Production-run leadtime Week 

FLPDTM Time of placing most recently de- A particular week 
livered production allocation rs 

FLQI Optimal ratio of EQQ to square root Unit 
of product demands 

FLR Reaction leadtime Week 

FLRI Review interval Week 

FLKIJK Variance correction factor for None 
demands and production leadtime 
to be applied to production lead= 
time demands 

FLZK Branch interim reorder leadtime; Week 
LSK + LRI/2 + LW/2 

FLS Leadtime, shipping, mean Week 

FLSK In=-transit leadtime, branch Days 

FMU oystem capacity constraint None 

FMUK Branch capacity constraint None 

FMDF Demand forecast nultiplier None 

FMBG Branch goals, multiplier None 

FMFG Plant goals, multiplier None 


Lop. 





FNOFR 


FRK K 
FZJK 


GENDEM 
HC 


ICON 
IC ONL 
ICOQN2 
ICON3 


ICONS 


ICON5 


IE OF 


Standard cost multiplier 
Leadtime multiplier, plant to 
warehouse 
Leadtime multiplier, warehouse to 
branch 
Number of production runs, floating 
point 
Freight rate 
Demands in branch interim order 
leadtime 
Subroutine to generate demands 
Holding cost due to rounding to 
pallets 
Parameter for RXCOMP subroutine 
Rule selector; 1-SRI, 2-other 
Not used 
Stockout policy selector, 
l- issue to zero balance 
2- do not make partial issues 
Round to cartons selector, 
l= do not round to full case lot 
2~- round to next higher case lot 
Round to pallets selector, 
l= do not 
2~ round if holding cost is less 
than handling costs 
Index 
End of file, fixed point mode 
Index 
Index for item 
Index used in subroutine predict 
Index used in subroutine predict 
Index for branch 


Index for month and production order 


Index and parameter for subroutine 
QUEULP 


Index parameter for subroutine QUEUEB 


Number of items 

Number of branches 

Order count, system 

Order count, product 

Order count, item 

Order count, item-location 

Order delivery date 

On-hand, item-location 

Shipping frequency 

On order, item-location, cumulative 
Order quantity 

Order count, branch, mean 

Page number 

Subroutine to place branch orders 
Plant goal 

Unit selling price, branch 
Production order allocation to item 


136 


None 
None 


None 
None 


Cents /Cwt 
Unit 


None 
Dollars 


None 
None 
None 
None 


None 


None 


None 
None 
None 
None 
None 
None 
None 
None 
None 


None 
None 
None 
Orders 
Orders 
Orders 
Orders 
Week 
Unit 
Week 
Unit 
Unit 
Orders 
None 
None 


Days (of Supply) 
Dollars/Unit 


Unit 


POAW J 


POBS J 
rObD 
Po 


PRIJK 


PREDICT 
QI 

Q1Z 
Q122 
QIPR 
QIPR 


QJ I 
QK I 


QORD 
QODD 


He 


QUEUVEB N 
QPOAS J,I 
QPOAW J,I 
QPOBS I 
i 
M 


9 
] 
QPODD 
QUEUE P 
QRIJK 

R 

RBI 


RGX 
RI 


RIJ 
RIJJ 
RIJK J,K 
RJK 

RJ1 

ALd 

RLTI 

RLTJ 


RLTJK 


Production order allocation to 
warehouse 

Production batch size 

Production order delivery date 

Picking cost due to roundg to 
cartons 

Probability of order between produc-=- 
tion allocations 

Subroutine to predict demands 

Economic prodcn order qty 

Computational factor for QI 

Minimum of alternate values of BOQ 

Computational factor for RIJK, XIJK 

The number of orders corresponding 
to the adjusted BOQ quantity 

Order item index in subroutine queue 

Order branch index in subroutine 
queue 

Order quantity in subroutine queue 

Order delivery date in subroutine 
queue 

Subroutine for queuing branch orders 

Production item allocatn in queue 

Production whse alloctn in queue 

Production batch size in queue 

Production order date in queue 

Subroutine for queuing production 
order 

Modified EOQQ value 

Scaling factor; POAS/ORD 

Requisition objective level, branch 
only 

Ratio for alternate trigger level 

Requisition objective level (goal), 
product 

Requisition objective level (goal), 
item 

Requisition objective level (goal), 
item 

Requisition objective level (goal), 
item-location 

Maximum of RZZJK and XIJK plus lead- 
time demands 

Renuisition objective level, ware- 
house only, item 

Ratio of losses to total demands, 
system 

Ratio of losses to total demands, 
product 

Ratio of losses to total demands, 
item 

Ratio of losses to total demands, 
item-location 
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Unit 


Product unit 
Week 
Dollars 


None 


None 
Preduct unit 
Froauct: univ 
Product unit 
Product unit 
None 


None 
None 


Unit 
A specific week 


None 

Unit 

Unit 

reouuct. whit 

A specific week 
None 


Product unit 
None 
Product unit 


None 
Product unit 


Unit 
Unit 
Unit 


Unit 


JUnat 


None 


None 


None 


None 


SOWJ 
SOWJ K 
SQ 
SQI 
SQP 


SQWI2 
SR 


ORP 


STI 
OTIK 


SYI 
SYJ 


SYJIK  J,K 


SOWK K 


Computational factor for requisition 
objective 


Stockturn 
Stockturn 
Stockturn 
Stockturn 
Stockturn 


ratio, 
rab. oO, 
ratio, 
ratio, 
ratio, 


Ssysten 
product 
item 
branch 


item=Location 


Run number 

Subroutine to compute SRI rules 

Alternate value of requisition 
objective 

Minimum of alternate values of RIJK 

Subroutine to read input data 

Size (item) code 

System demands 

System demds less warehouse demds 

System demds in leadtime 

Sys demds in leadtime, branches only 

System goal all products, and 
locations 

System goal, less warehouse 

System goal, warehouse 

Inventory position, product 

Inventory position, item 

Inventory position, item-location 

Product unit conversion factor 

Stockout count, system 

Stockout count, product 

Stockout count, item 

Stockout count, item-location 

Cumulative ord qty 

Stockout weeks, system 

Stockout weeks, product 

Stockout weeks, item 

Stockout weeks, item-lLocation 

Economic order quantity, system 

Economic order quantity, product 

Economic order quantity, branches 
only 

Square of Wilson's EQQ, product 

System requisition objective, all 
prods + locs 

System requisition objective, 
branches only 

Coverage factor, product 

Coverage factor, location 

Average inventory position, 

Average inventory position, 

Average inventory position, 

Average inventory position, 
location 

Stockout weeks, branch 

Time index 


system 
product 
item 
item- 
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Unit 


None 
None 
None 
None 
None 
None 
None 
Unit 


Unit 
None 
None 
Product 
Product 
Product 
Product 
Product 


Product 
Proauct 
Product 
Una 
Unit 
Product 
None 
None 
None 
None 
Product 
Week 
Week 
Week 
Week 
Preduct 
Product 
Product 


Product 
Product 


Product 


Week 
Week 
Product 
Product 
Unit 
Unit 


Week 
None 


unit 
UnLY 
unit 
unit 
unit 


unit 


unit 
unit 


unit/Unit 


nat 


Unaet 
unit 
unit 


unit 
unit 


unit 


unit 
unit 





UL 
ULI 
ULJ 
ULJK 
UPCJ 
VARDF 
VDIJ1 
VDFR 
VDIJK 
VDS 
VDSI 
VDSJ 
VDSK 
VDSJ K 


VLP 
VL2K 
VSIT 
VSYI 
VS YJ 
VSYK 
VSYJK 
VesiIT 
VZD 
VabI 
VZDJ 
VZDK 
VLIK 
VDZJK 


V2 
V3 


V4 
V5 
v8 
v2 


Coverage provided by QI (BAR) plus 
one half of reaction time 

Temporary storage location for OHJK 

Temporary storage location for DSIJK 

Shelf life 

Interval between time of placing 
most recently delivered production 
allocation and present time 

Units loss, system 

Units loss, product 

Units loss, item 

Units loss, item-location 

Units per carton 

Function to compute variance 

Rel-variance of demands, whs 

Probability of interim order 

Rel-variance of demand, item-location 

Value of demand satisfied, system 

Value of demand satisfied, product 

Value of demand satisfied, item 

Value of demand satisfied, branch 

Value of demand satisfied, item- 
location 

Rel-variance of prodcn leadtime 

Rel-variance of branch leadtime 

Value of inventory in-transit 

Value of ave inventory, product 

Value of ave inventory, item 

Value of ave inventory, branch 

Value of ave inventory, item-location 

Value of ave inventory, system 

Value of ave inventory in-transit 

Value of total demands, system 

Value of total demands, product 

Value of total demands, item 

Value of total demands, branch 

Rel-Variance of branch leadtime 

Rel-variance of demands, item- 
location 

Control knob - system requisition 
level in weeks of demand except 
for central whse under SRI rules 

Control knob - system requisition 
level in weeks of demand includ- 
ing central whse under SRI rules 

Control knob - difference between 
V2 and Wl 

Control knob = difference between 
V3 and W2 

Control knob = expected maximum on 
hand level, all branches 

Control knob = expected average on 
hand level, all branches 
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Week 


None 
None 
Month 
Week 


Product unit 
Product unit 


Unit 
Units 


Units/Carton 


None 
None 
None 
None 
Dollars 
Dollars 
Dollars 
Dollars 
Dollars 


None 

None 

Dollars 
Dollars 
Dollars 
Dollars 
Dollars 
Dollars 
Dollars 
Dollars 
Dollars 
Dollars 
Dollars 
None 

None 


None 


None 


None 
None 
None 


None 





V10 


Vill 


WPSOT 
WFSOJ 
WPSOJK 
W1 


W2 


XCI J 


0 
XDIJK J,K 


a 

XIJ J 
mak Oy K 
XJ1 J 
XJK 

Mork J 


XRIJK 
XRIJZ J 


XX1J J 


XXIJZ J 


ZDIJ 
ZDIJK 
ZVC 
ZVCK 


af 


Control knob - same as V8 except 
warehouse is excluded 

Control knob = same as V9 except 
warehouse is excluded 

Temp storage location 

Unit weight 

Weeks per stockout, 

Weeks per stockout, 

Weeks per stockout, item 

Weeks per stockout, item-location 

System requisition level in weeks 
of demand except for central 
whse under alternate rules 

System requisition level in weeks 
of demand including central 
whse under alternate rules 

Temporary storage location 

Dealer const (unit plus shipping 
plus receiving) 

Demands over shelf life, item 


system 
product 


Demands over shelf life, item-loca- 


tion 
Reorder trigger product 


item 


level, 


Reorder trigger 


Reorder 
Reorder 
Reorder 

other 


floor 


mand 


trigger 
trigger 
trigger 


than warehouse 
Interim order leadtime 


level, 
level, 
level, 
level, 


item-location 
warehouse 
item-location 


times reorder 


multiplier (FK2) 

Requisition objective level ata 
branch in weeks of item location 
demands 

Requisition objective level at ware- 
house in weeks of total item de- 


Reorder level for item in weeks of 
total item demand 
Reorder level for warehouse in weeks 


of total item demand 


Temporary storage location 
Total demands item 
Total demands item-location 
Total of CH, OC, CL; system 
Total of CHK, CBOK, CLK; branch 
Obsolescence cost (not used in 
present form of simulator) 
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None 
None 


None 
Pounds /Unit 
Week 
Week 
Week 
Week 
Week 


Week 


None 
Dollars/Unit 


Unit 
Unit 


Product unit 
Unit 
Upac 
Unit 
Unit 


None 


Week 


Week 


Week 
Week 


None 

Unit 

Unit 
Dollars 
Dollars 
Cents/Unit 
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